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CHAPTER  1 

INTRODUCTION 

Amplification  by  a  Modulated  Electron  Beam 

The  ability  of  a  modulated  electron  beam  to  propagate  space-charge 
waves  was  first  investigated  by  Hahn  [I]  and  Ramo  [2].   Passing  a  beam  of 
electrons  through  a  gap,  across  which  an  r-f  electric  field  is  applied, 
causes  the  electrons  to  be  decelerated  or  accelerated  depending  on  the 
phase  of  the  electric  field  at  the  time  of  passage.   The  accelerated 
electrons  tend  to  catch  up  with  the  electrons  which  left  the  gap  at  an 
earlier  time,  but  were  decelerated  by  the  applied  electric  field.   This 
grouping  of  electrons  is  referred  to  as  bunching.   It  is  assumed  that 
the  charge  density  of  an  unmodulated  electron  beam  is  zero,  as  there 
are  equal  numbers  of  electrons  and  positive  ions.   The  positive  ions 
result  from  the  ionization  of  the  residual  gases  by  the  electron  beam. 
In  the  modulated  electron  beam,  the  charge  density  in  the  region  where 
the  electrons  are  bunched  is  negative;  whereas,  in  the  region  between 
bunches,  the  charge  density  is  positive.   Because  of  this  difference 
in  charge  density,  the  space-charge  forces  between  the  electrons  move 
electrons  from  the  region  of  high  density  into  the  region  of  low  density. 
Since  regions  of  lower  electron  density  appear  both  ahead  of  and  behind 
a  region  of  high  density,  some  electrons  are  accelerated  with  respect  to 
the  average  beam  velocity  while  others  are  decelerated.   To  a  stationary 

1 


observer,  there  would  appear  to  be  two  waves  propagating  with  velocities 
slightly  higher  and  lower  than  the  average  beam  velocity.   As  the  bunches 
of  electrons  spread  apart  in  the  axial  direction  under  the  influence  of 
space-charge  repulsion,  new  bunches  of  electrons  form  from  the  decel- 
erated electrons  of  one  bunch,  and  the  accelerated  electrons  of  the 
following  bunch  formed  by  the  r-f  field  a  short  time  later.   A  second 
gap  is  positioned  at  some  point  along  the  beam  where  a  new  bunch  of 
electrons  has  formed.   This  large  grouping  or  bunching  of  electrons  means 
a  sizeable  component  of  r-f  current  is  present  in  the  beam.   This  current 
will  induce  a  similar  component  of  r-f  current  in  the  output  cavity,  pro- 
vided the  cavity  is  tuned  to  the  r-f  frequency  at  which  the  beam  is  ex- 
cited.  The  cavity  may  be  tuned  instead  to  a  harmonic  of  the  fundamental 
frequency,  in  which  case  the  device  is  considered  a  frequency  multiplier 
rather  than  an  amplifier.   The  resultant  voltage  produced  in  the  output 
cavity  is  an  amplified  version  of  the  input  excitation.   The  analyses  of 
Hahn  and  Ramo  shows  that  the  optimum  position  for  the  output  gap  is  one- 
fourth  of  a  plasma  wavelength,  where  the  plasma  frequency  is  the  fre- 
quency at  which  the  bunched  electrons  oscillate  under  the  effect  of 
space-charge  forces.   The  bunching  process  takes  place  without  further 
Influence  by  external  fields  if  the  beam  Is  enclosed  in  the  region  be- 
tween cavities  by  a  conducting  cylinder  called  a  drift  tube.   The  two 
cavities  with  the  drift  tube  separating  them  are  the  essential  elements 
of  the  microwave  device  called  the  klystron.   Because  of  the  use  of 
tuned  resonant  cavities,  the  klystron  has  an  Inherent  bandwidth  limita- 
tion, and  another  device  called  the  t ravel i ng-wave  tube  has  been 


developed.   Essentially  the  drift  tube  is  replaced  by  a  helix,  which 

supports  an  electromagnetic  wave  which  propagates  with  an  axial  veloc- 
ity less  than  the  speed  of  light.   Thus   the  electromagnetic  wave  can 

interact  with  the  space-charge  wave  on  the  beam  in  a  continuous  fashion. 

This  helix  eliminates  the  need  for  the  tuned  output  cavity,  which  has  a 

limited  bandwidth,  since  the  amplified  signal  is  taken  from  the  helix 
di  rectly. 

Small-Signal  Analysis 

The  action  of  the  electron  beam  may  be  described  by  specifying 
how  the  velocity  v(z,t)  and  current  i(z,t)  modulations  fluctuate.   A 
first  approximation  for  the  solutions  of  velocity  and  current  may  be 
obtained  by  a  process  referred  to  as  small-signal  analysis.   Essentially 
the  equations  are  derived  by  utilizing  the  equations  for  the  space-charge 
force  on  an  electron;  the  continuity  of  current;  the  relations  between 
velocity,  current,  and  charge  density;  and  Maxwell's  equations.   In 
this  approach  it  is  assumed  that  the  time  variation  is  sinusoidal,  and 
the  equations  are  linearized  by  discarding  the  product  of  a-c  terms. 
The  separation  of  a-c  and  d-c  terms  then  yields  equations  which  are 
functions  of  position  (z)  alone,  since  the  variation  with  time  (eJU)t) 
cancels  in  the  a-c  equation. 

Space-Charge-Wave  Amplifiers 

Tien  and  Field  [3]  reexamined  the  principles  of  klystron  operation, 
and  found  that  further  amplification  could  be  obtained  if  the  drift  tube 


was  divided  into  sections  of  appropriate  length,  each  having  a  d-c 
potential  different  from  the  adjacent  sections.   Additional  energy  is 
available  for  increasing  the  amplitude  of  velocity  modulation  over  that 
which  would  occur  in  a  single  section  of  drift  tube.   This  device, 
called  the  velocity-jump  amplifier,  was  the  first  of  a  series  of 
devices  classified  as  space-charge-wave  amplifiers.   Others  in  this 
class  are  the  rippled-wall  amplifier,  where  the  sections  of  the  drift 
tube  take  on  different  diameters  without  changing  potentials,  and  the 
r ippled-stream  amplifier,  where  the  drift  tube  is  unchanged,  but  the 
beam  diameter  fluctuates.   The  former  device  was  investigated  by 
Birdsall  and  Whinnery  [k] ,    and  the  latter  by  Birdsall  [5]  and  Mihran  [6]. 
A  general  analysis  of  a  space-charge-wave  amplifier  which  incorporates 
all  these  possible  alternatives  has  been  performed  by  Peter,  Bloom  and 
Ruetz  [7]  using  small-signal  theory.   It  was  found  that  maximum  ampli- 
fication resulted  if  the  drift-tube  potential  and  the  drift-tube  diameter 
were  increased  while  the  beam  cross  section  was  decreased  at  the  point 
where  the  a-c  velocity  was  minimum.   At  a  velocity  maximum  all  para- 
meters were  then  returned  to  their  original  values.   One  interesting 
approach  for  investigating  the  operation  of  space-charge-wave  amplifiers 
is  a  construction  by  Bloom  and  Peter  [8]  of  a  transmission-line  analog 
which  has  a  standing-wave  pattern  that  represents  the  amplitude  pattern 
of  the  space-charge  wave.   The  approach  is  based  on  the  small -signal 
analysis  which  predicts  that  the  maximum  current  and  velocity  amplitudes 
interchange  every  quarter  plasma  wavelength.   A  general  discussion  of 
the  linearized  small-signal  analysis  as  applied  to  the  various  microwave 


amplifiers  is  presented  in  books  by  Beck  [9],  Harman  [10],  and  Hutter  [ll]. 
The  small-signal  theory  is  presented  in  Chapter  2,  and  the  nonlinear  analy- 
sis Is  the  subject  of  the  following  chapters. 

Finite  Beam  Considerations 

Since  beams  are  of  finite  diameter  and  located  close  to  metallic 
conductors,  i.e.,  the  walls  of  the  drift  tube,  it  is  apparent  that  the 
space-charge  forces,  which  cause  the  bunching  and  debunching  of  the  beam, 
will  be  reduced  because  of  the  fringe  effect  of  the  electric  field.   This 
fringe  effect  increases  as  the  beam  radius  (b)  approaches  that  of  the 
drift  tube  (a),  and  is  more  prominent  in  beams  having  a  small  normalized 
beam  radius  (^eb) .   The  effect  of  reducing  the  space-charge  force  is  a 
reduction  in  the  plasma  frequency  (cj  ),  i.e.,  the  frequency  at  which  the 
beam  bunching  takes  place.   This  modification  is  obtained  by  multiplying 
the  plasma  frequency  by  a  parameter  called  the  plasma  reduction  factor 
(R.).   Hahn  made  the  first  analysis  of  finite  beams  in  his  original 
report,  and  discovered  that  an  infinite  number  of  modes  of  space-charge- 
wave  propagation  can  exist.   Ramo  suggested  that  for  a  reasonable  approxi- 
mation only  the  fundamental  mode  need  be  considered.   More  rigorous  sum- 
mation techniques  demonstrate  that  the  effect  of  the  higher  order  modes 
Is  most  pronounced  in  beams  which  effectively  fill  the  drift  tube  and 
have  high  Seb  values.   Branch  and  Mihran  [12]  have  solved  the  trans- 
cendental equations  which  result  from  the  Hahn-Ramo  theory  for  the  radial 
propagation  constant  needed  to  determine  the  reduction  factor.   Curves 


are  presented  with  Rj  plotted  against  either  Seb  or  b/a,  with  the 
remaining  variable  as  a  parameter,  for  several  beam  geometries.   The 
plasma  reduction  factor  affects  the  wavelengths  of  the  standing-wave 
envelopes  of  the  velocity  and  current  modulations,  and  the  amplitude 
of  the  current  is  also  inversely  proportional  to  the  reduction  factor. 
The  plasma  reduction  factor  is  derived  under  the  assumptions  of  small- 
signal  theory,  and  one  of  the  problems  involved  in  large-signal  analysis 
has  been  the  applicability  of  the  reduction  factor  to  large-signal 
situations.   Mlhran  [13]  has  demonstrated  experimentally  that  replace- 
ment of  cj  by  R,w   in  the  extrapolated  small-signal  theory  leads  to 
errors  in  both  magnitude  and  location  of  current  maxima. 

Extensions  to  Small -Signal  Theory 

The  small-signal  theory  correctly  predicts  the  fundamental  opera- 
tion of  these  microwave  devices  for  small  ratios  of  buncher  voltage  (V  ) 

m 

to  beam  voltage  (Vn) .   As  the  signal  level  is  raised,  experimentally 
observed  phenomena  such  as  harmonic  generation,  saturation  of  the  fund- 
amental current,  phase  shift,  and  growth  of  the  second  harmonic  component 
are  not  even  predicted,  much  less  analyzed,  by  the  small-signal  approach. 
Several  approaches  have  been  conceived  to  extend  the  linear  analysis  into 
the  nonlinear  region  of  operation.   These  schemes  are  discussed  in  the 
following  sections. 

Extension  of  Webster's  theory 

One  of  the  earliest  analyses  of  klystron  operation  was  performed 


by  Webster  []k].      The  Webster  solution  for  the  fundamental  component  of 
current  is 

|f  =  2I0J,(X0).  (1.1) 

where 

X0  =  ballistic  bunching  parameter  |l(a/2)Sez, 

|Q  =  d-c  beam  current, 

J]  =  Bessel  function  of  1st  order  and  1st  kind, 

|l  =  gap  coupling  coefficient, 

a     =   voltage  modulation  index  (Vm/V0), 

V  =  amplitude  of  signal  voltage  applied  to  buncher, 

Vq  =  d-c  beam  voltage, 

Se  =  u/v0. 

This  solution  is  based  on  the  assumption  of  no  space-charge  interaction, 
and  hence  is  a  ballistic  analysis.  The  Webster  theory  has  been  modified 
to  account  for  space-charge  effects  by  comparing  (1.1)  with  the  solution 
derived  by  the  small -signal  space-charge  wave  theory, 

|f  =  l0(a/2)nBe2  sin(Bpz)/(Bpz),  (1.2) 

where  B  =  w  /vQ.   The  solutions  given  by  (1.1)  and  (1.2)  yield  identical 
results  for  small  values  of  a.   The  Bessel  function  may  be  replaced  by 
one-half  of  its  argument  for  small  values  of  the  argument.   Calling 
the  argument  X,  and  comparing  (1.1)  and  (1.2)  under  the  assumption 
that  X  is  smal 1  gives 


X  =  H(ar/2)Bez  s in(Bpz)/(SpZ) .  (1.3) 

This  solution  for  the  bunching  parameter  when  space-charge  effects  are 
accounted  for  is  substituted  into  (1.1),  and  used  as  the  solution  for  lr 
regardless  of  the  value  of  X.   There  is  no  mathematical  justification 
for  this  substitution,  and  it  is  used  simply  because  an  extrapolated 
solution  is  better  than  no  solution  at  all.   Comparison  of  the  ballistic 
bunching  parameter  (XQ)  defined  in  (1.1)  with  the  solution  obtained  in 
(1.3)  shows  that  the  bunching  parameter  including  space-charge  effects 
is  related  to  XQ  by  the  relation 

X  =  XQ  sin(Bpz)/(Spz).  (1.4) 

The  term  sIn(B  z)/(B  z)  is  called  the  debunching  factor.   As  a    increases 
the  modified  Webster  theory  indicates  that  the  current  does  not  increase 
proportionally,  i.e.,  the  fundamental  component  of  current  saturates. 
Although  the  modified  solution  has  been  extrapolated  past  the  point  of 
electron  crossover,  the  predicted  ballistic  paths  of  the  electrons  are 
incorrect.   The  modified  theory  fails  to  account  for  a  change  in  the 
direction  of  the  space-charge  force  past  crossover,  and  for  a  change  in 
the  magnitude  of  the  force.   Moreover,  the  theory  Is  one-dimensional, 
i.e.,  the  fringing  effect  of  the  finite  beam  Is  neglected.   Mihran  [13] 
traces  the  development  of  a  further  modification  to  Webster's  theory 
for  predicting  the  electron  paths  past  crossover.   This  theory,  developed 
by  Roe,  corrects  the  first  two  errors  of  the  extrapolated  theory;  however, 


the  analysis  Is  still  one-dimensional.   In  order  to  approximate  the  effect 
of  the  fringe  field,  the  plasma  reduction  factor  is  introduced  into  the 
equations  by  replacing  B  by  RjB  .    Mihran  demonstrates  that  klystron 
design  based  on  this  approximation  results  in  significantly  less  than 
optimum  efficiency.   He  cites  a  typical  large-signal  case  where  the 
efficiency  is  20  per  cent  less  than  an  empirical  optimum. 

Electron  disk  approach 

Tien,  Walker,  and  Wolontis  [15]  devised  the  idea  of  replacing  the 
electron  beam  by  a  set  of  charged  disks  which  have  radii  equal  to  the 
beam  radius,  and  are  evenly  distributed  along  the  axis  of  the  beam. 
They  applied  this  electron  disk  approach  to  the  traveling-wave  tube  for 
small  values  of  the  gain  parameter  (C) ,  and  later  Tien  [16]  extended 
the  analysis  to  larger  values  of  C.   Webber  [17,  18]  applied  the  disk 
model  to  the  klystron,  and  Rowe  [19]  applied  it  to  the  travel ing -wave 
tube  using  a  space-charge  field  based  on  an  electron  distribution  in 
time  rather  than  in  space.   Although  the  disk  approach  includes  the 
effect  of  space  charge  in  an  approximate  manner,  the  idea  of  replacing 
a  plasma  of  electrons  by  16  or  2k   charged  disks  is  questionable. 
Webber's  klystron  computations  do  not  fit  Mihran's  [13]  experimentally 
derived  values  of  "peel -off  angle  anywhere  over  a  wide  range  of  Beb 
values.   The  "peel-off1  angle  is  that  normalized  distance  of  the  drift 
tube  for  which  the  normalized  r-f  current  is  one-half  db  below  its 
ballistic  maximum.   This  model  does,  however,  predict  the  growth  of 
the  second  harmonic  which  Mihran  [20]  verified  experimentally. 
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Paschke's  successive  approximation 

An  alternative  approach  has  been  applied  by  Paschke  [21],  who  con- 
sidered the  electron  beam  as  a  fluid  just  as  was  done  for  the  linear 
case  by  Hahn  and  Ramo.   Maxwell 's  equations  are  applied  to  the  model, 
and  the  resulting  nonlinear  partial  differential  equation  is  solved 
without  discarding  products  of  a-c  terms  or  assuming  e-*   variation  in 
time  as  the  small-signal  analysis  does.   The  technique  for  solving  the 
nonlinear  equation  is  referred  to  as  the  method  of  successive  approxi- 
mations.  It  is  assumed  that  the  dependent  variable  may  be  expressed 
as  a  power  series  In  the  voltage  modulation  coefficient  (a).   The  space- 
charge-wave  equation  is  then  grouped  into  terms  according  to  degree. 
The  first-order  solution,  which  is  proportional  to  a,    is  obtained  by 

retaining  only  the  first-degree  terms.   This  solution  is  then  used  to 

2 
obtain  the  second-order  solution  which  is  proportional  to  cc    ,  and  is 

found  by  retaining  terms  In  the  general  equation  which  are  proportional 
to  cc    .      This  process  Is  repeated  until  it  is  assumed  that  the  additional 
terms  have  a  negligible  contribution,  i.e.,  the  power  series  is  con- 
vergent.  In  the  approach  used  by  Paschke,  a  third-order  theory  was 
considered  to  be  sufficient  to  investigate  the  phenomena  of  saturation, 
phase  shift,  and  second  harmonic  generation.   Thus  the  solution  is  of 

the  form 

3 

.-£l„.  (1-5) 


where  in  is  proportional  to  a    .   The  approach  has  a  fundamental  assumption, 
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which  is  shared  with  the  disk  model,  that  the  electric  field  does  not 
vary  over  the  cross  section  of  the  beam.   It  has  been  shown  by  Olving  [22] 
and  Mihran  [13]  that  this  assumption  is  valid  for  thin  beams  (Seb  <  1). 
The  main  limitation  to  the  theory  Is  the  fact  that  it  is  only  valid  for 
velocities  which  are  single-valued  functions  of  space,  i.e.,  no  electron 
overtaking  is  permitted.   In  succeeding  reports,  Paschke  considered  finite 
beams  with  fringe  fields  in  order  to  evaluate  the  second  harmonic  com- 
ponent [23],  saturation  of  the  fundamental  and  third  harmonic  generation 
[24],  and  phase  shift  of  the  fundamental  [25].   The  phase-shift  problem 
has  also  been  investigated  by  Engler  [26]  without  using  some  of  the 
approximations  utilized  by  Paschke  [25].   Olving  [22,  27]  applied  the 
method  of  successive  approximations  to  a  velocity-modulated  electron 
beam  and  obtained  an  infinite  series  solution.   The  results  included 
the  effect  of  higher  order  modes,  whereas  Paschke  had  assumed  that  only 
the  fundamental  mode  was  generated.   Olving  demonstrated  that  Paschke 's 
intuitive  Idea  of  using  plasma  reduction  factors  in  the  nonlinear  region 
Is  justified  for  thin  beams  where  higher  order  modes  of  the  harmonics 
may  be  neglected.   Since  this  restriction  also  applies  to  the  approxi- 
mation of  constant  electric  field  over  the  beam  cross  section,  the  method 
of  successive  approximations  introduces  no  new  restrictions  other  than 
the  one  of  no  electron  overtaking.   Although  this  method  does  not  con- 
sider overtaking,  Paschke  [2k]    showed  that  his  method  fits  the  experi- 
mental values  of  "peel-off"  angle  defined  by  Mihran  [13]  very  well  in 
the  region  where  B  b  <.   1,  i.e.,  in  the  region  where  the  model  is  valid. 
With  the  method  of  successive  approximations  It  Is  possible  to 
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investigate  right  up  to  the  point  of  electron  overtaking  with  some 
assurance  that  the  results  can  be  verified  by  experiments. 

Scope  of  the  Present  Investigation 

Paschke's  method  of  analysis  is  applied  to  the  space-charge  waves 
propagating  in  a  space-charge-wave  amplifier  which  combines  the  features 
of  the  velocity-jump  and  rippled-wall  amplifiers.    Beam  diameter  changes 
are  not  considered  as  it  is  relatively  difficult  to  produce  this  vari- 
ation physically,  whereas  the  other  two  changes  can  be  done  with  com- 
parative ease.   Some  of  the  phenomena  investigated  in  the  first  region 
include  the  saturation  of  the  fundamental  current  as  the  bunching  volt- 
age is  increased,  phase  shift,  second  harmonic  generation,  and  the  shift 
in  locations  of  the  velocity  and  current  maxima  and  minima  from  the  pre- 
dicted small-signal  locations.   Investigation  in  the  second  region 
includes  the  effect  of  different  ratios  of  drift-tube  potentials,  the 
effect  of  changing  the  drift-tube  radius  between  sections,  and  the  location 
of  the  transition  between  the  two  regions  in  order  to  cause  maximum  current 
or  velocity  variations  for  fixed  sets  of  beam  potentials  and  drift-tube 
geometries.   These  changes  will  be  investigated  with  the  thought  in  mind 
of  how  the  current  excited  in  an  output  cavity  can  best  be  increased  for 
some  fixed  value  of  a   by  properly  utilizing  the  velocity-jump  and 


'rlclsaac  [28]  applied  this  nonlinear  analysis  to  a  device  with  an 
alternate  series  of  gaps  and  drift  tubes;  however,  no  provision  was  made 
for  applying  d-c  potentials  across  the  gaps,  or  modifying  the  drift-tube 
diameters. 
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r ippl ed-wal 1  principles  In  combination.   The  proper  lengths  of  the  drift 
tube  sections  in  order  to  optimize  the  output  current  will  also  be  ex- 
amined, as  will  the  saturation  effect  which  limits  the  output  current 
as  or  Increases. 

All  symbols  utilized  in  this  paper  are  defined  in  Appendix  A,  and 
all  constants  defined  in  the  derivation  of  the  space-charge-wave  equation 
are  evaluated  in  Appendix  B. 


CHAPTER  2 

SMALL-SIGNAL  ANALYSIS  OF  A  SPACE-CHARGE-WAVE  AMPLIFIER 

Synopsis 

The  smal 1 -signal  solutions  of  the  Hahn-Ramo  theory  are  applied  to 
a  space-charge-wave  amplifier  composed  of  three  drift-tube  sections. 
The  first  and  the  third  sections  have  identical  properties,  while  the 
middle  section  possesses  both  a  different  potential  and  a  different 
drift-tube  radius.   Investigation  includes  both  increasing  and  decreas- 
ing the  potential  and  radius  of  the  second  drift  tube,  relative  to  the 
values  of  these  parameters  in  the  first  region.   The  gain  is  found  to 
be  proportional  to  the  ratio  of  the  plasma  reduction  factors  and  the 
three-fourths  power  of  the  voltage  ratio.   Maximum  gain  results  if  the 
lengths  of  the  drift  tubes  are  some  multiple  of  one-quarter  of  a  re- 
duced plasma  wavelength. 

Small-Signal  Equations 

The  small -signal  equations  developed  by  Hahn  and  Ramo,  as  given  by 
Hutter  [11],  are 

j  (ut  -  B  z) 
i  -  [i.  cos(6qZ)  -  j(to/wq)  (iQ/vo)  v.  sin(SqZ)]  e 

(2.1) 
14 
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v  =  [v.  cos(Sqz)  -  j(uq/u>)  (vQ/io)  I,  s!n(Bqz)]  eJ((*  "  6eZ) 


(2.2) 


where  the  \_   subscript  indicates  an  Initial  value  at  the  entrance  to  the 
region  described  by  the  equations,  and  j  =v-l.   The  other  variables 
were  defined  in  connection  with  (1.1)  and  (1.2),  and  are  tabulated  for 
easy  reference  in  Appendix  A.   These  equations  are  derived  on  the 
assumptions  that  the  variation  with  time  is  given  by  e^   ,  that  a-c 
terms  are   small  compared  with  d-c  quantities,  and  that  the  theory,  which 
is  derived  for  a  beam  of  infinite  cross  section,  can  be  applied  to  beams 
of  finite  diameter  by  use  of  the  plasma  reduction  factor  as  tabulated  by 
Branch  and  Mihran  [12].   A  plasma  reduction  factor  curve  for  a  beam  of  ra- 
dius (b)  inside  a  drift  tube  of  radius  (a)  is  shown  in  Fig.  2.1.   A  list  of 
the  properties  required  of  the  plasma,  in  order  to  derive  the  small- 
signal  equations,  is  stated  in  Chapter  3  in  connection  with  the  deri- 
vation of  the  nonlinear  space-charge-wave  equations.   Equations  (2.1) 
and  (2.2)  are  valid  for  any  region  in  which  the  d-c  velocity  (v  )  is 
constant,  such  as  inside  a  section  of  drift  tube. 


Initial  Conditions 

The  velocity-jump  amplifier,  developed  by  Tien  and  Field  [3], 
utilizes  a  change  in  d-c  potential  between  drift-tube  sections  to 
further  amplify  the  current  and  velocity  modulation.   In  order  to 
use  (2.1)  and  (2.2)  to  describe  the  amplification  process,  it  is 
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2.4 


B.b 


NORMALIZED  BEAM  RADIUS  (B  b) 


Fig.  2.1  -  Space-charge  reduction  factors  for  cylindrical  beams 
In  cylindrical  drift  tubes. 
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necessary  to  determine  what  happens  to  the  magnitude  of  the  a-c  current 
and  velocity  as  the  beam  passes  through  the  transition  point.   Hutter 
develops  a  proof,  using  the  Lorentz  force  equation,  which  demonstrates 
that  the  a-c  current  is  continuous,  whereas,  the  a-c  velocity  undergoes 
a  transition  proportional  to  the  ratio  of  the  d-c  velocities  in  the  two 
regions. 

Ibdt)  "  ia(zt)  (2-3) 

vb<zt>  =  <Vvb>*  va<zt>'  (2-*> 

where  the  subscripts  a_  and  b_  refer  to  the  first  and  second  regions 
respectively,  and  zt  is  the  position  at  which  the  transition  is  located. 
A  model  of  the  amplifier  is  shown  in  Fig.  2.2. 

Before  the  beam  enters  the  input  cavity,  located  at  z0,  there  is 
no  a-c  excitation.   A  signal  voltage  V  cos(wt)  velocity  modulates  the 
beam  at  the  input  cavity.   The  initial  values  of  a-c  current  and 
velocity  are  derived  in  Chapter  3  in  connection  with  the  nonlinear 
analysis.   The  small-signal  solutions  are 


l,(z0)  -  0  (2.5a) 


Vi<2o>  "  <a/2)  Voa'  (2'5b) 


where  a    is  the  ratio  of  signal  voltage  (V  )  to  d-c  beam  voltage  (V  ). 
The  equations  for  the  first  drift-tube  region  are  found  by  substituting 
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(2.5)  into  (2.1)  and  (2.2). 

ia=  -  Ha/2)    (U/Uqa)  iQ  sin(Sqaz)  (2.6) 

va  "  (a/2>  voa  cos(Bqaz>  <2'7) 

The  exponential  term,  which  carries  only  phase  shift  information,  has 
been  omitted  from  (2.6)  and  (2.7),  and  is  also  omitted  in  the  rest  of 
the  equations  in  this  chapter.   To  derive  any  instantaneous  solution, 
it  is  necessary  to  insert  the  exponential  term,  and  take  the  real  por- 
tion of  the  product.   Only  envelope  magnitudes  are  considered,  therefore 
this  exponential  term  may  be  omitted. 

Small -Signal  Analysis 

If  the  first  drift  tube  is  one-fourth  of  a  reduced  plasma  wavelength 
long  (SQaz  =  90°),  then  the  magnitude  of  the  a-c  variations  at  za  are 

!.<*.)  =  (0/2)  (u»/Uqa)  i0  (2.8a) 

va(za)  =  0.  (2.8b) 

Equations  (2.8)  are  utilized  to  obtain  entrance  conditions  for  the 
second  drift-tube  region.   Solutions  for  the  second  region  are  obtained 
by  substituting  (2.8)  into  (2.1)  and  (2.2). 


20 


ib  =    (a/2)    (co/u     )    iQ  cos(B     z)  (2.9a) 


j    (a/2)    (to  ,/u)     )    v        sin(B     z)  (2.9b) 

J  v    qb     qa        ob  qb  v  ' 


If   the    length  of   the   second   drift    tube    is   B  bzb  =  90°,    then 

ib(zb)    =  0  (2.10a) 

vb(zb)    =    (a/2)    (u)qb/uqa)    vQb.  (2.10b) 

The  entrance  conditions  for  the  third  section  are   obtained  by  comparing 
(2.10)  and  the  requirement  given  by  (2.*+). 

i,(zb)  =  0  (2.11a) 

Vzb>  =  <VVa>J  Vb(2b)  <2',,b> 

The  solutions  for  the  third  region  are   obtained  by  substituting  (2.11) 
into  (2.1)  and  (2.2).   Since  the  third  region  has  the  same  properties 
as  the  first  region,  the  parameters  use  the  a  subscript. 

ic  =  -  j  (a/2)  («/«qa)  («qb/Wqa)  (Vb/Va)  i0  *in(6qaz)   (2- 12a) 
vc  =    (a/2)  <VV  ^b/Va)"  vob  cos(Bqaz)      (2J2b) 

The  ratio  of  the  maximum  a-c  velocity  in  the  third  region  (vcm)  to  the 
maximum  a-c  velocity  in  the  first  region  (vam)  is 
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v   /v   =  (V./V  )*  (R./R  )  (u  ,/u  )  (v  ,/v   ).  (2.13) 

cm  am   v  b  a'   x  b  a'  x  pb  pa'  x  ob  oa'  v     ' 

The  d-c  velocity  is  proportional  to  (V0)  ,  and  the  plasma  frequency  is 

inversely  proportional  to  the  square  root  of  the  d-c  velocity.   Thus 
(2. 13)  reduces  to 

v™/vam  "  (Rh/RJ  (VVJ3A-  (2-|Z+) 
cm  am     b  a    b  a 


The  ratio  of  the  current  maxima  is  identical  to  (2.14). 

Examination  of  (2.14)  indicates  that  the  best  amplification  is  ob- 
tained if  the  values  of  the  potential  and  plasma  reduction  factor  in  the 
second  drift-tube  region  are  increased.   The  value  of  Ru  can  be  increased 
by  Increasing  the  drift-tube  radius,  as  shown  in  Fig.  2.1.   If  the  size 
of  the  drift  tube  is  not  increased,  the  value  of  R^  is  less  than  R 
because  the  normalized  beam  radius  (3  .b)  decreases  as  the  drift-tube 
potential  is  increased.   The  variations  in  the  envelopes  of  the  current 
and  velocity  modulations  are  shown  in  Fig.  2.3,  for  a  choice  of  drift 
tube  and  beam  parameter  values  dealt  with  on  a  nonlinear  basis  in  later 
chapters. 

The  nonlinear  investigation  includes  only  the  first  two  regions. 

The  value  of  maximum  velocity  in  Region  C,  which  occurs  just  after  the 

second  transition  (z^) ,  is  obtained  by  multiplying  the  value  of  the 

velocity  just  before  the  second  transition  by  the  ratio  of  the  d-c 

velocities  (v  k/v  ).   This  yields  adequate  information  for  demonstrat- 
v  ob  oa         '         n 

ing  the  velocity  saturation  effect.   The  modified  Webster  theory, 
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discussed  in  Chapter  1,  can  be  utilized  to  predict  the  a-c  current  present 

in  the  output  cavity  located  a  quarter  plasma  wavelength  beyond  the  second 

transition.   The  current  prediction  is  obtained  by  replacing  a  in  (1.3) 

by  a  new  a',  where  a      is  equal  to  2v c_/v   ,  and  substituting  (1.3)  into 

(1.1)  in  place  of  X0.   In  addition  Bp  is  replaced  by  R]Bp  to  account  for 

the  effect  of  the  fringe  field. 

An  expression  similar  to  (2.14)  results  If  the  first  transition  is 

made  at  one-half  of  a  plasma  wavelength  (B   z  =  180°).   The  entrance 
K  3   v  qa  a 

conditions  for  the  second  drift  region  are 

l,(za)  =  0  (2.15a) 


v;(za)  =  (a/2)  (VV^  voa.  (2.15b) 


The  new  solutions  for  the  second  drift  space  are  found  by  applying  (2.15) 
to  (2.1)  and  (2.2). 

ib  =  -  j  (a/2)  (w/u,qb)  (Va/Vb)  i0  sin(Bqbz)     (2.16a) 

vb=  (a/2)  (Va/Vb)2  voa  cos(B  z)  (2.16b) 


The  second  drift  tube  must  be  one-fourth  of  a  reduced  plasma  wavelength 
long  in  order  to  allow  the  second  transition  to  occur  at  a  current 
maximum.   The  entrance  conditions  for  the  third  region  are 

i.(zb)  =  (a/2)  (u)/coqb)  (Va/Vb)  lQ  (2.17a) 

v,(zb)  -  0.  (2.17b) 
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The  solutions  in  the  third  region,  assuming  once  again  that  this  region 
possesses  the  same  properties  as  the  first  region,  are  obtained  by  sub- 
stituting (2.17)  into  (2.1)  and  (2.2). 


ic  =  (a/2)  (ai/u^)  (Va/Vb)  i0  cos(6qaz)      (2.18a) 

v  =  -  j  (a/2)  (R  /R.)  (V  /VU)3A  v   sin(S  z)   (2.18b) 
c  a   b    a   b      oa      qa 


The  variations  in  the  envelopes  of  the  current  and  velocity  modulations, 
with  the  first  transition  occurring  at  the  point  of  maximum  velocity 
variation,  are  shown  in  Fig.  2.*+.   The  ratio  of  the  maximum  current  or 
velocity  In  the  third  region  to  the  corresponding  maximum  value  in  the 
first  region  is 

i rm/iam  =  v  /v   =  (Ra/Rj  (V„/V,)3A.  (2.19) 

cm  am    cm  am   v  a   b'  v  a  b'  *    ' 

It  is  apparent  from  examining  (2.19)  that  in  order  to  obtain  maximum 
amplification  it  is  necessary  to  decrease  the  drift-tube  potential  and 
radius  in  the  second  region.   This  is  just  the  opposite  approach  from 
that  of  the  first  analysis;  however,  comparison  of  (2.14)  and  (2.19) 
indicates  that  the  degree  of  amplification  is  the  same.   The  only 
difference  comes  in  the  length  of  the  drift-tube  sections.   With  the 
first  transition  occurring  at  a  current  maximum,  all  drift  tubes  are 
one-fourth  of  a  reduced  plasma  wavelength  long;  whereas,  in  the  second 
approach,  the  first  and  third  sections  must  be  one-half  of  a  plasma 
wavelength  long  in  order  to  obtain  the  desired  current  amplification. 
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Conclusions  Based  on  Small -Signal  Theory 

The  small-signal  theory  predicts  that  the  amplification  is  inde- 
pendent of  the  degree  of  voltage  modulation,  I.e.,  it  is  independent 
of  a.      The  experimental  fact  that  amplifiers  actually  reach  saturation 
quite  quickly  as  a    increases  cannot  be  ascertained  from  this  small- 
signal  theory.   This  theory  also  predicts  that  the  best  amplification 
is  obtained  if  the  voltage  is  decreased  at  a  velocity  maximum,  and  in- 
creased at  a  velocity  minimum.   Because  the  small-signal  solution  for 
a-c  velocity  has  a  minimum  value  of  zero,  no  reduction  in  the  a-c 
velocity  occurs  if  the  beam  voltage  is  increased  at  the  point  of  zero 
a-c  velocity.   If  the  transition  is  located  at  some  other  point,  then 
the  amplitude  of  the  a-c  velocity  is  reduced.   These  simple  results 
are  contrasted  in  Chapter  6  with  the  results  obtained  by  a  nonlinear 
theory.   Among  other  things,  it  is  shown  that  the  best  amplification  is 
not  obtained  by  placing  the  d-c  velocity  transitions  at  positions 
corresponding  to  multiples  of  a  quarter  plasma  wavelength,  and  that 
saturation  effects  quickly  come  into  play  as  the  signal  amplitude  is 
increased. 


CHAPTER  3 

DERIVATION  OF  SPACE-CHARGE-WAVE  EQUATIONS  USING 
THE  METHOD  OF  SUCCESSIVE  APPROXIMATIONS 

Synops  i  s 

Following  the  approach  used  by  Paschke,  the  electron  beam  is  con- 
sidered as  a  plasma,  where  the  effect  of  the  individual  electron  is 
ignored.   The  nonlinear  partial  differential  equation  which  describes 
the  polarization  of  an  electron  is  derived  by  examining  the  current 
flow  in  an  infinitesimal  beam  element.   Linear  equations  are  obtained 
by  the  method  of  successive  approximations. 

Assumpt  ions 

In  order  to  simplify  the  analysis  the  following  assumptions  are 
made : 

1.  The  electron  velocity  is  small  compared  with  the  velocity  of 
1  ight. 

2.  The  r-f  electric  field  component  in  the  axial  direction 
is  approximately  constant  over  the  beam  cross  section. 

3.  The  effects  of  temperature,  particle  collisions,  and  potential 
depression  caused  by  space  charge  are  igncred. 

k.       Ion  mass  is  assumed  to  be  infinite,  hence  the  ions  are  not 
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affected  by  r-f  fields. 

5.  The  d-c  charge  density  of  the  ions  (lp0l)  cancels  the  d-c  charge 
density  of  the  electrons  (p0) . 

6.  The  velocity,  in  keeping  with  the  polarization  coordinate  system, 
is  a  single-valued  function  of  space,  i.e.,  electron  overtaking 
does  not  occur. 

7.  The  effect  of  the  transverse  electric  field  associated  with  the 
plasma  surrounded  by  a  conducting  cylinder  is  accounted  for  by 
using  the  appropriate  plasma  reduction  factor. 

8.  All  variables  are  assumed  to  be  functions  only  of  the  independent 
variables  time  (t)  and  axial  position  (z). 

9.  An  infinite  magnetic  field  confines  the  motion  of  the  electron 
beam  to  the  axial  direction. 

Derivation  of  the  Space-Charge-Wave  Equations 

Several  relations  between  position,  current,  velocity,  and  charge 
density  are  needed  in  order  to  derive  the  nonlinear  space-charge-wave 
equations  which  describe  the  electron  plasma.   The  charge  density  is 
given  by 

pt  =  po  +  p(z.t).  (3.1) 

The  convection-current  density  is  related  to  the  charge  density  and 
velocity  by 
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!t  =  'o  +  '  =  (Po  +  P)(vo  +  v)'  (3.2) 

where 

io=Povo-  (3-3^ 

The  a-c  current  density  and  charge  density  are  related  by  the  principle 
of  continuity  of  current, 

-  di/dz  =  dp/dt.  (3.4) 

Applying  Newton's  notation  to  partial  differential  equations  gives 

df/dz  a  f',  (3.5a) 

and 

df/dt  =  f,  (3.5b) 

where  £  is  any  function  of  axial  position  (z)  and  time  (t).   Thus  applying 
(3-5)  to  (3.4)  gives 

-  i'  =  p.  (3.4) 

The  velocity  at  any  point  in  the  plasma  is  a  function  of  the  effect  of 
the  space-charge  forces  at  that  point.   Using  Newton's  law  of  motion  this 
relationship  is  given  by 

v  +  (vQ  +  v)  v'  =  -  tjE.  (3.6) 
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There  is  an  associated  electric  field  transverse  to  the  z-axis  when' 

ever  an  electron  beam  of  finite  size  is  in  the  vicinity  of  a  metallic 

drift  tube.   This  transverse  electric  field,  which  is  time  varying, 

indicates  the  presence  of  a  displacement  current  (M)  per  unit  length 

of  beam.   This  coupling  current  enters  into  the  continuity  equation  for 

the  total  current  flowing  through  an  infinitesimal  length  of  beam  as 

shown  in  Fig.  3- 1 • 

A  M  dz 


Fig.  3-1  -  Current  flow  in  an  infinitesimal  beam  element 


Summation  of  the  currents  shown  in  Fig.  3-1  gives 


i   +  eo(E)   +  M/A  =  0. 


(3.7) 


At  this  point  a  new  variable  (z  )  is  introduced.   This  variable  repre- 
sents the  polarization,  or  displacement,  of  an  electron  from  its  position 
in  the  unmodulated  beam.   This  choice  of  variable  is  in  keeping  with  the 
polarization  system  of  coordinates  established  by  Bobroff  [29].   The 
polarization  is  related  to  the  electric  field  due  to  space-charge  by 
the  relation 
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-  Po  zp  =  €o  E  +  f,  (3.8) 

where  f_  is  the  coupling  flux  which  represents  the  presence  of  the  trans- 
verse electric  field.   The  coupling  flux  is  related  to  the  coupling 
current  (M)  by  the  relation 

(f)'  =  M/A.  (3.9) 

Substitution  of  (3-8)  and  (3-9)  into  (3-7)  gives 

i'  -  P0(ip)'.  (3-10) 

Integration  of  (3-10)  with  respect  to  z  yields 

1  "  Po  V  (3.11) 

Comparison  of  (3.10)  and  (3.*0  in  terms  of  i   gives 


P=  -PQ(ip)  •  (3.12) 

Integration  of  (3-12)  with  respect  to  time  yields 


P  =  "  P0zp  •  (3-13) 

The  velocity  is  obtained  by  substituting  relations  for  i  and  p  into  (3.2) 


v  =  V°  'P  +  fr  (3.14) 

zp 


32 


An  equation  for  velocity  is  obtained  by  solving  (3-8)  for  E  and  substi- 
tuting the  result  into  (3-6). 


v  +  (vQ  +  v)  v  =  H(pozp  +  f)/€o  (3.15) 


Substituting  (3-1^)  into  (3-15)  results  in  an  equation  involving  only 
the  polarization  (zp)  of  the  electron.   The  final  form  is  a  nonlinear, 
second-order  partial  differential  equation  of  fourth-degree. 


zp"  +  2(zp)'/v0  +'zp/v02  -  il(p0zp  +  f)/e0v02 

+  2[  2p(zp)'  -  Zp'zp  -  (Zp)'zp'zp  ]/v02 

+  [z  (z  ')2  +  z  "  (z  )2]/v  2  +  2[z  "z  -  (z  ) '  z  ']/v 
Pv  P       P    P     °       P  P     P    P   o 

+n(P0zp  +  f)  t3zp'  -  3(zp')2  +  (Zp')3  ]/£0v02  =  0     (3.16) 


The  velocity  is  related  to  the  polarization  (z  )  through  (3.1*0i  and 
the  current  density  is  related  to  the  polarization  through  (3.11). 


Method  of  Successive  Approximations 

The  solution  to  (3.16)  is  obtained  by  the  technique  known  as  the 
method  of  successive  approximations.   In  this  approach  it  is  assumed 
that  the  polarization  may  be  expressed  as 


(3.17) 
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where  it  is  assumed  that  z   is  proportional  to  an .      Three  equations  are 
formulated  by  substituting  (3-17)  into  (3.16)  and  separating  terms  by 
powers  of  a.   The  variables  i,  y_,  and  f_  are  also  related  to  a   by  series 
simi lar  to  (3-1 7) - 

Substitution  of  (3.17)  into  (3.16)  leads  to  a  first-order  equation 
of  the  form 

(z,)"  +  2(z,)7v0  +  2,/v02  -  (Ve0)(P0z,  ♦  fi)/v02  «  0.(3.18) 

Equation  (3.18)  has  a  solution  (zj)  which  is  directly  proportional  to  a. 
This  solution  should  be  identical  to  the  Hahn-Ramo  small-signal  solution, 
which  is  also  a  linear  function  of  a.      The  Hahn-Ramo  solution  is  obtained 
from  (3.18)  provided, 


Po*l  +  fl  ~  Rl2Pozl-  (3.19a) 


This  is  another  way  of  saying  that  the  small-signal  plasma  reduction 
factors  may  be  utilized  to  account  for  the  effect  of  the  fringe  field. 
Following  the  technique  used  to  obtain  (3.19a),  Paschke  has  approximated 
the  second-order  and  third-order  fringe  effects  by  the  relations 

Poz2  +  f2  =  <  R0^2  >   Po*2  <3J9b> 

Poz3  +  f3  =  <  R,^  >  poz3.  (3.19c) 

The  notation  <  Rq  2  >  ^ol lows  the  examPle  of  Engler,  and  indicates  that 
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R2,  the  value  of  reduction  factor  associated  with  2B  b,  Is  utilized  when 
seeking  the  second  harmonic  portion  of  the  second-order  solution,  whereas 
Rq  =  0  is  used  with  the  d-c  portion  of  the  second-order  solution.   The 
third-order  solution  has  a  component  of  fundamental  frequency  and  a 
third  harmonic  portion.   Hence,  Ri  is  used  with  the  fundamental  frequency 
portion,  and  R,  with  the  third  harmonic.   The  value  of  R,  is  determined 
by  examining  the  reduction  factor  curve  of  fixed  b/a  at  the  point  3Beb. 
Olving  [22]  has  shown  that  the  application  of  plasma  reduction  factors 
is  valid  for  thin  beams  (Beb  <    1). 

Substituting  (3.19a)  into  (3.18)  gives 


z,"  +  2  (*,)'/vo  +  kyv/  +  (R,Bp)2  z,  =  0,  (3-20) 

where  the  plasma  wave-number  (B  )  is  defined  as 


Bp  =  Up/Vo  =,7niPol/e0/v0.  (3.21) 


Three  linear  partial  differential  equations  are  obtained  from  (3.16)  by 
substituting  (3.17)  and  (3-19)  into  (3-16)  and  separating  according  to 
powers  of  a. 

z,"  +  2(z,)7v0  +'z,/v02  +  (BpR,)2  Z]  =  0  (3.22a) 

z2"  +  2(z2)'/vo  +-z2/vo2  +  Bp2  <  R022  >  z2  =  3(RlSp)2  z]Z]' 
-  2  [z,  (i,)'  -  ZjZ,']^2  -  2  [2]"  z,  -  (z,)'  Z|']/vo        (3.22b) 
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z3  +  2(z3)  /vo  +  z3/v/  +  Bp«  <  R,^  >  z3  = 

2  [  Z]'z2  +  Zj'z,  -  z,(z2)'  -  i^i,)']/^2 

+  [  2(Z])'  z/z,  -  Zjtz,')2  "  z/'tz,)2]/^2 

+  2  [(z,)'  z2'  +  (z2)'  z,'  -  z,"  z2  -  z2"  z,]/v0 

+  V  C3«l  '  <  R022>  22  "  3^  R,2>3  >  z1(z]')2  +  3z2'  <  R,23  >  z,] 

(3.22c) 

The  velocity  may  also  be  expressed  as  a  series  solution 

3^ 
v  -  £  vn,  (3.23) 

n=l 

where  vn  is  proportional  to  ccn . 

Expanding  the  denominator  of  (3-1*+)  into  a  power  series  yields 

v  -  (v0  zp'  +  zp)  [1  +  zp'  +  (zp')2  +  (zp')3  +  ...  ].      (3.24) 

Equations  for  the  first-,  second-,  and  third-order  velocity  terms  are 
obtained  by  substituting  (3-17)  into  (3-24)  and  comparing  coefficients 
of  powers  of  a. 

v,  =  vo  z,'  +  z,  (3.25a) 

v2  =  v,  z,'  +  z2  +  vQ  z2'  (3.25b) 

v,  =  v  z  '  +  z„  +  v.  z  '  +  v  z.  (3.25c) 

3    o  3    3    12    2 
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The  solution  for  the  electron  polarization  (z  )  is  obtained  by  first 
solving  the  linear,  homogeneous,  partial  differential  equation  (3.22a)  for 
the  first-order  solution  (zj).   This  solution  is  substituted  into  the 
right  hand  side  of  (3.22b)  in  order  to  evaluate  the  forcing  function.   The 
second-order  polarization  (z2)  is  obtained  by  solving  the  linear,  nonhomo- 
genous   equation  (3.22b).   The  third-order  polarization  is  obtained  by 
solving  (3.22c),  after  first  substituting  the  solutions  for  Z|  and  z2 
into  the  right  hand  side  of  (3.22c)  to  evaluate  the  forcing  function.   For 
this  investigation  a  third-order  approximation  is  considered  to  be  suffi- 
cient to  determine  the  saturation  effect.   Consequently,  the  electron 
polarization  (zp)  is  approximated  by  the  sum  of  Z| ,  z2,  and  z^ . 


CHAPTER  k 

GENERAL  SOLUTIONS  FOR  THE  SPACE-CHARGE-WAVE  EQUATIONS 

Synops  i  s 

In  this  chapter  the  general  solutions  for  the  first-,  second-,  and 
third-order  components  of  the  polarization,  current,  and  velocity  are 
derived.   Unlike  the  analyses  of  Paschke  and  Engler,  boundary  condi- 
tions are  not  substituted  upon  completion  of  each  order  of  the  solu- 
tion.  Consequently,  the  general  solutions  are  considerably  more  complex 
than  the  solutions  obtained  in  their  reports.   Both  Engler  and  Paschke 
restricted  their  analysis  to  the  investigation  of  a  single   drift  region, 
corresponding  to  Region  A  of  this  report.   Since  the  primary  investigation 
in  this  report  is  the  effect  of  drift-tube  potential  and  reduction  factor 
changes  on  the  amplification  process,  at  least  two  unique  drift  regions 
are  needed.   The  additional  labor  required  to  obtain  a  general  solution 
is  necessary  since  the  particular  solution  in  Region  B  is  obtained  by 
comparison  of  the  general  solution  obtained  in  this  chapter  with  the 
particular  solution  for  Region  A  obtained  in  the  next  chapter. 

First-Order  Solutions 

Since  the  small-signal  solution  to  the  space-charge-wave  equation 
has  a  time  dependence  given  by   eJ '    "  eZ,      a  solution  of  the  form 
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z,  =  F,(z)  cos(wt  -  Bez)  +  G,(z)  sin(ut  -  Bez)  (4.1) 

Is  chosen.   Equations  for  the  functions  Fj(z)  and  Gi(z)  are  obtained  by 
substituting  (4.1)  into  (3.22a). 


[F,"(z)  +  (R,B  )2  F,(z)]  cos(u)t  -  Sez) 


+  [G,  (z)  +  (R,B  )Z  G^z)]  sin(iot  -  B  z)  =  0       (4.2) 


Since  (4.2)  must  hold  for  all  values  of  time, 


F,"(z)  +  (R,Bp)2  F,(z)  =  0  (4.3a) 


G,  (z)  +  (R,6p)2  G,(z)  =  0.  (4.3b) 


The  product  R]B   is  called  the  reduced  plasma  wave  number  and  is  defined 
by  the  symbol  (B  ,).   The  solutions  to  (4.3)  are 


F,(z)  =  A  cos(Bq)z)  +  B  sin(Bq]z)  (4.4a) 


G,(z)  =  C  cos(Bq]Z)  +  D  sin(Bq)z).  (4.4b) 


The  general  solution  for  the  first-order  polarization  (zj)  is  obtained 
by  substituting  (4.4)  into  (4.1).   The  polarization  is  normalized  by 
multiplying  (4. 1)  by  B Q]a.   The  normalized  first-order  polarization  is 
given  by 
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h    ■  Bqla  zl  =  &   cos(Bq,2)  +  B  sin(Sq,z)]  cos(ut  -  B&z) 

+    [C  cos(6q)z)  +  5  sin(Bq]z)]  sin(u>t  -  Bez)  ,  (4.5) 

where  the  bar  (~)  indicates  that  the  parameter  has  been  multiplied  by 

B  ,  .   The  subscript  a  of  B„.   indicates  that  this  is  the  value  of  B_i 
q  I  a  r   —     q  i  a  q  i 

in  Region  A.   The  absence  of  a  second  subscript  for  the  other  terms  in 
(4.5)  simply  indicates  that  the  terms  are  applicable  for  either  region. 
The  particular  solutions  for  Region  A  or  Region  B  have  double-subscripted 
terms . 

The  first-order  current  (i,)  is  found  by  substituting  (4.5)  into 
(3.11). 

i,/i   =  BE  [C  cos(B  ,z)  +  5  sin(B  ,z)]  cos(cot  -  B  z) 
lo  v  q  I  x  q  I  '      N      e 

-  BE  [A  cos(Bqlz)  +  B  sin(Sqlz)]  sin(ut  -  Bez) , 


(4.6) 


where 


BE  =  VBqla' 


The  first-order  velocity  is  found  by  substituting  (4.5)  into  (3.25a). 

v]/v0  =  BQ.1  [B  cos(Bq]z)  -  A  sin(Bq)z)]  cos(ut  -  Bez) 

+  BQ1  [6  cos(Bq|z)  -  C  sin(Bq]z)]  sin(u)t  -  Bez) ,   (4.7) 

where 

BQJ  =  B  ,/B  .  . 
^     ql   qla 
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Second-Order  Solutions 


The  substitution  of  (4.5)  into  (3.22b)  determines  the  forcing 
function  for  the  second-order  polarization  equation. 

zZ"   +  2  (z2)  7v0  +  z2/v02  +  Bp2  <  rq22  >  z%   =  ^         (Zf>8) 

where 

ff2  =  K,  cos(2Bq]z)  +  K2  sin(2Bq]Z) 

+  [K3  +  K^  cos(2B  ,z)  +  K5  sin(2Bq]z)]  cos  2(ut  -  Bez) 

+  [K6  +  K  cos(2B  ]Z)  +  Kg  sin(2B  z)  ]  sin  2(ujt  -  B  z) 

The  constants  K]  through  Kg  are  defined  in  Appendix  B.   The  general 
solution  for  z2  is  of  the  form 

z2  =   F0^    +  ?2^    cos   2(CJt    "  Bez)    +  Mz)    sin   2(£Jt    *  Bez)'    ^-,0) 

The  reduction  factor  Rq  is  associated  with  the  d-c  portion  of  (4.8) ,   and 
R2  with  the  second  harmonic  portion.   Three  equations  are  obtained  by 
substituting  (4.10)  into  (4.8)  and  comparing  coefficients  of  the  same 
function  of  time. 

Fo"(z)  +  <RoV2  F0<z)  =  Kl  cos(2Bq]z)  +  K2  sin(2Bqlz)      (4.11a) 
F2"(z)  +  (R26p)2  F2(z)  =  K3  +  K^  cos(2Bq]z)  +  K^   sin(2Bq]Z)  (4.11b) 
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G2"(z)  +  (R2Sp)2  G2(z)  =  K6  +  K?  cos(2Bq]z)  +  Kg  sin(2Bq]z)  (4.11c) 

Since  Rq  =  0,  the  solution  to  (4.11a)  is  found  by  integrating  twice  with 
respect  to  z. 

FQ(z)  =  -  [K,  cos(2Bq)z)  +  K2  s in(2Bq| z)  ]/(2Bq , ) 2  +  E  z  +  H, 

C+.12) 

where  E  and  H  are  constants  of  integration.   The  complementary  function 
for  (4.11b)  is 

F2(z)c  =  I  sin(Bq2z)  +  J  cos(Bq2z),  (4.13) 

where 

Bq2  =  R2Bp.  (4.14) 

The  form  of  the  particular  integral  is  determined  by  the  forcing  function 
port  ion  of  (4. lib). 

F2(z)p  =  C,  +  C2  cos(2Bq,z)  +  C3  sin(2Bqlz)       (4.15) 

The  arbitrary  constants  Cj ,  C2 ,  and  Co  are  evaluated  by  substituting 
(4.15)  into  (4.11b)  and  comparing  coefficients. 


F2(2)p  =  [K3  +  g  (K4  cos(2Bq]Z)  +  «5  s i n(2Bq ,z) ) ]/(Bq2)2,    (4.16) 


where 


g  =  (i  -  *♦  s2)-\ 
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and  s  =  Bql/Bq2. 


The  general  solution  to  (4.11b)  Is  the  sum  of  (4.13)  and  (4.16). 
Comparison  of  (4.11b)  and  (4,11c)  indicates  that  the  form  of  the 
solution  to  (4.11c)  is  identical  to  (4.16)  provided  K^  is  replaced 
by  Kgfetc.   The  complementary  function  involves  two  new  constants  of 
integration  called  L  and  M.   The  general  solution  for  z.  is  the  sum 
of  the  general  solutions  for  FQ(z) ,  F2(z),  and  G2(z).   As  in  the  case 
of  the  first-order  solutions,  the  second-order  solutions  are  normalized. 

z"2  =  E  z  +  H  -  SqIa  [K,  cos(2Bq]z)  +  K2  s  in(26q]  z)  ]/(2Bq,  )2 

+  [T  sin(Bq2z)  +  J  cos(6q2z)  +  Bq]a  {K3  +  g  Kk   cos(2Bq]z) 
+  g  K5  sin(2Bq,z))/(Bq2)2]  cos  2(ut  -  Bez)  +  [C  sln(Bq2z) 
+  M  cos(Bq2z)  +  Bq|a  {Kg  +  g  K?  cos(2Gq]z)  +  g  Kg  sin  (2Bq,z)} 
/(Bq2)2]  sin  2(ut  -  Bez)  (4.17) 

The  second-order  current  is  found  by  applying  (4.17)  to  (3.11). 

I2/IQ  =  2  BE  [M  cos(Bq2z)  +  L  sin(Bq2z)  +  Bq]a  (Kg  +  g  K? 

•cos(2Bq]Z)  +  g  Kg  s ?n(2Bq] z)  )/(Bq2)2]  cos  2(wt  -  Bez) 

-  2  BE  [J  cos(Bq2z)  +  I  sin(Bq2z)  +  Bq)a  (K3  +  g  K4 

•cos(2B  ,z)  +  g  K5  s in(2B  . z) }/(B  2)2 ]  sin  2(ut  -  Bez) 

H  M       H  (4.18) 
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The  second-order  velocity  Is  found  by  substituting  (4.5) ,    (4.7),  and 
(4.17)  Into  (3.25b). 

v2/vQ  =  E  +  BE  BQ1  P]0/2  +  (BQJ)2  P2/4  -  [K2/(6Bq])]  cos(2Bqlz) 
+  [K,/(6Bq])]  sln(2Bq]z)  +  [BQ2  T  cos(Bq2z)  -  BQ2  J 
•sin(Bq2z)  -  Kg/Bg  +  (1  +  2s2)  g  Uj  cos(2Bqlz)  -  Kk 
•sin(2Bq]z))/(3Bq,)]  cos  2(wt  -  Bez)  +  [BQ2  L  cos(Bq2z) 

-  BQ.2  R  sin(Bq2z)  +  K^/B&   +  (1  +  2s2)  g  {K3  cos(2Bq]Z) 

-  K?  sin  (2Bq]z)}/(3Bq])]  sin  2(wt  -  Bez)  (4.19) 


Third-Order  Solutions 

The  forcing  function  for  the  third-order  polarization  equation  is 
evaluated  by  substituting  solutions  for  zj  and  z2  Into  (3.22c). 

23"  +  2(z3)'/v0  +  z3/v02  +  <  R,23  >  Bp2  z3  =  [Kg  cos(Bq]Z) 
+  K,0  cos(b]Bq2z)  +  K,,  cos(b2Sq2z)  +  K]2  cos(3Bqlz)  +  K,3  sin(Bq]z) 
+  K]k   sin(b|Bq2z)  +  K]5  sin(b2Bq2z)  +  K|6  sin(3B  ]2)]  cos(u)t  -  ^z) 
+[K,7  cos  (Bq]z)  +  K]8  cos(b]Bq2z)  +  K,g  cos(b2Bq2z)  +  K20 
•cos(3Bqlz)  +  K2]  sin(Bq]z)  +  «22  sin(b]Bq2z)  +  K23  sin(b2Bq2z) 
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+  K2k   sin(3Sq]z)]  sin(ut  -  Sez)  (4.20) 

Only  that  portion  of  the  forcing  function  which  is  of  fundamental 
frequency  has  been  retained.   The  other  portion,  which  is  a  third 
harmonic  function  of  time,  and  hence  is  associated  with  the  generation 
of  a  third  harmonic  component  of  polarization,  is  omitted.   The  reason 
for  this  omission  is  simply  that  the  problem  is  concerned  with  the 
fundamental  frequency,  and  harmonic  components  are  carried  only  when 

they  are  required  to  evaluate  solutions  which  involve  fundamental 

2  2 

frequency  components.   The  notation  <  R.  ,  >  is  replaced  by  Rj   because 

of  the  omission  of  the  third  harmonic  component. 

The  form  of  the  general  solution  for  the  third-order  polarization 

is 

z3  =  F3(z)  cos(wt  -  Bez)  +  G3(z)  s?n(u)t  -  Bez) .        (4. 21 ) 

The  equations  needed  to  determine  F,(z)  and  G3(z)  are   obtained  by  sub- 
stituting (4.21)  into  (4.20)  and  comparing  coefficients  of  the  same 
function  of  time. 


F3"(z)  +  (R,Bp)2  F3(z)  =  A,  (4.22a) 


G3"(z)  +  (R,Bp)2  G3(z)  =  B,,  (4.22b) 


where  A]  represents  the  coefficient  of  cos(wt  -  Bgz)  and  B]  of  sin(ujt  -  Bez) 
as  given  in  (4.20).   The  solutions  of  the  complementary  functions  of  (4.22) 
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F3(z)c  =  P  cos(Bq]z)  +  0_  sin(Bq]Z)  (4.23a) 

G3(z)c  =  R  cos(Bq]z)  +  S  sin(Bq]z).  (4.23b) 

The  solutions  for  the  particular  integrals  are  determined  by  the  forcing 
function  in  (4.20).  Since  cos(B  .z)  is  part  of  both  the  forcing  function 
and  the  complementary  function,  the  general  form  of  the  particular  inte- 
gral must  be  modified  to  include  the  term  Ci  z  cos(B  .z).  A  similar 
substitution  is  made  for  the  sin(B  ,z)  term.  The  various  coefficients 
are  evaluated  by  substituting  the  general  form  into  (4.20)  and  comparing 
with  the  forcing  function.   The  general  solutions  for  F3(z)  and  G-j(z)  are 

F3(z)  =  [P  -  K]3z/(2Sq,)]  cos(Bq|z)  +  [Q  +  K9z/(2BqI)]  sin(Bq]z) 
+  £K,0  cos(b,Bq2z)  +  K]k   sin(b]Bq2z)]/(b3Bq22) 

'  [K11  COS^b2Bq22)  +K15Sin(b2V)]/(b4Bq22) 

-  [K)2  cos(3B  z)  +  K]6  sin(3B  z)]/(8B  2)        (4.24) 


G3(z)  =  [R  -  K2]z/(2Bq|)]  cos(Bq]z)  +  [S  +  K]7z/(2B  ,)]  sin(Bq]z) 
♦  [K]8  cos(b]Bq2z)  -K22sin(bIBq2z)]/(b3Bq22) 

-  [K)9  cos(b2Bq2z)  +  K23  sin(b2Bq2z)]/(b4Bq22) 

-  [K20cos(3Bqlz)  +  K2k   sin(3Bq)z)]/(86ql2).         (4.25) 
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The  general  solution  for  the  normalized  third-order  polarization  is 

z3  =  F3(z)  cos(ut  -  Bez)  +  G3(z)  sln(ut  -  tSQz) ,  (4.26) 

where  the  bar  (")  indicates  that  the  term  is  multiplied  by  Gq]a. 

The  third-order  current  is  obtained  by  substituting  (4.26)  into 
(3.11). 

i3/i0  -  BE  G3(z)  cos(wt  -  Sez)  -  BE  F3(z)  s»n(wt  -  BQz)  (4.27) 

The  third-order  velocity  is  determined  by  substituting  zj,  z2,  Zo  ,  Vi, 
and  V2  into  (3.25c) . 

VVo  =  ^3<2)  +  Fa(z)  +  Fc(z)]  cos(t0t  "  Sez) 

+  [G3'(z)  +  Fb(z)  +  Fd(z)]  sin(uit  -  Bqz)  ,     (4.28) 

-  '        _  i 
where  F3  (z)  and  G3  (z)  are  the  derivatives  with  respect  to  z  of  F,(z) 

and  G3(z)  respectively,  and  the  functions  Fa(z),  F|j(z),  Fc(z),  and 

F  .(z)  are  defined  in  Appendix  B. 

In  this  chapter  the  general  solutions  to  the  space-charge-wave 

equations  are  determined.   Boundary  conditions  are  applied  in  the  next 

chapter  at  the  input  gap  to  determine  the  particular  solutions  in 

Region  A.   In  the  following  chapter,  boundary  conditions  are  applied 

at  the  point  where  the  drift  tube  properties  change  to  determine  the 

particular  solutions  in  Region  B.   The  solutions  and  variables  which 
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apply  in  one  particular  region  will  be  denoted  by  a  second  subscript. 
Thus  the  first-order  polarization  in  Region  A  wi 1 1  be  denoted  as  Zj  , 
while  in  Region  B,  the  solution  is  denoted  as  z,.. 


CHAPTER  5 

PARTICULAR  SOLUTIONS  FOR  REGION  OF  FIRST  DRIFT  TUBE 

Synops  is 

In  this  chapter  particular  solutions  are  obtained  from  the  general 
solutions  of  Chapter  k   by  applying  boundary  conditions  at  the  input 
cavity  located  at  z  =  0.   The  solutions  for  this  initial  drift-tube 
region,  known  as  Region  A,  are  utilized  as  entrance  conditions  for 
obtaining  solutions  for  Region  B.   The  information  obtained  from  the 
solutions  in  Region  A  pertains  to  a  basic  two-cavity  klystron  amplifier, 
driven  at  low  signal  levels,  with  the  nonlinear  beam  effects  taken  into 
account.   The  phenomena  investigated  include  current  saturation,  shift 
in  the  locations  of  velocity  and  current  maxima  and  minima  as  signal 
voltage  increases,  phase  shift,  and  second  harmonic  current  generation. 
All  constants  identified  with  Region  A  wi  1  1  bear  an  a_  subscript. 
Numerical  data  used  to  plot  the  various  curves  presented  in  this  chapter 
were  obtained  by  use  of  a  Fortran  compiler  applied  to  the  IBM  709  computer. 

Boundary  Conditions 

A  signal  voltage  (Vs)  Is  applied  to  a  pair  of  closely-spaced  grids 
located  at  z  =  0.   The  signal  voltage  is  given  by 
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kS 


Vs(t)  =  Vm  sin(u)t).  (5.1) 

where  w  is  the  angular  frequency  at  which  the  voltage  varies.   A  previously 
unmodulated  beam  of  electrons  passes  through  the  pair  of  grids.   The  effect 
of  this  a-c  voltage  on  the  electron  beam  is  the  establishment  of  velocity 
modulation.   The  extent  of  the  effect  is  determined  by  applying  the 
principles  of  continuity  of  current  and  conservation  of  energy.   The  usual 
klystron  assumptions  [30]  discussed  in  the  following  sections  are  applied 
at  the  input  gap. 

Continuity  of  current 

The  principle  of  continuity  of  current  is  used  to  determine  the 
initial  current  at  z  =  0.   Since  the  beam  is  unmodulated  before  entering 
the  gridded  gap,  the  current  modulation  upon  exit  is  assumed  to  be  neg- 
1 igible.   Hence, 

Tj (0)  =  l2(0)  =  i3(0)  =  0.  (5.2) 

This  conclusion  is  based  on  the  assumption  that  the  two  grids  are  so 
closely  spaced  that  there  is  no  time  for  electron  bunching  to  occur, 
i.e., the  transit  time  is  zero.   Using  the  assumption  of  closely-spaced 
grids,  it  may  be  assumed  that  the  beam  coupling  coefficient  (|i)  is  unity. 
For  gaps  with  finite  grid  spacing  the  coefficient  is  always  less  than 
one.   Since  it  is  assumed  that  the  transit  time  is  zero,  it  follows  that 
the  initial  polarization  (z_),  or  displacement  from  the  undisturbed  posi- 
tion, is  also  zero.   This  condition  is  used  to  complete  the  second-order 
polarization  solution. 
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Conservation  of  energy 

The  velocity  of  the  electron  beam  upon  exit  from  the  input  cavity  is 
obtained  by  applying  the  principle  of  conservation  of  energy  at  i  =  0, 
along  with  the  assumption  that  p.  =  1  . 


Transposing, 


S  ince 


(m/2)  (va2  -  voa2)  =  -e  ^sinCut)  (5-3) 


va2  =  voa2  +  2  ^  Vm  sln(wt).  (5.4) 


voa2  eMVa,  (5.5) 


where  Va  is  the  potential  of  the  drift  tube  in  Region  A,  (5.4)  may  be 
written  as 

va2  =  2  n   Cva  +  Vm  sln(ut)].  (5.6) 

The  velocity  in  Region  A  is  obtained  by  factoring  Va,  taking  the  square 
root,  and  replacing  n/2  i)  Va  by  v   . 


va  =  voasM  +  (Vm/Va)  sin(uit)  (5-7) 

The  ratio  Vm/Va  is  by  definition  the  voltage  modulation  index  (a). 

a=Vm/Va.  (5-8) 
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Application  of  the  binomial  expansion  to  (5-7)  yields 

va  =  voa  [1  +  (a/2)  sfn(u)t)  -  (a/2)2  sin2(wt)/2  +  (a/2)3 

•  sin3(urt)/2  +  ...  ].  (5.9) 

Completing  the  square  and  cube  in  (5-9)  gives 


va  "  voa  ^  +  (a/2)  sin(ojt)  -  {a/k)2   {1  -  cos(2ut)} 

+  (a/4)3  (3  sin(ut)  -  sin(3ut))  +  ...  ].  (5-10) 


By  applying  Paschke's  method  of  successive  approximations  the  beam 
velocity  may  be  written  as 

va  ■  voa  +  vla  +  v2a  +  v3a>  (5-") 

where  vna  is  proportional  to  an.   The  various  velocity  terms  may  be 
evaluated  at  z  =  0  by  comparing  (5.  10)  and  (5-11)  by  powers  of  a. 

v,a(0)  =  (a/2)  sin(ut)  (5.12a) 

2 
v2a(0)  -  -  (a/k)      [1  -  cos(2wt)]  (5.12b) 

v3a(0)  =  (a/4)3  [3  sin(wt)  -  sln(3ut)],  (5.12c) 

where 

7na(°)  3  vna(°)/voa-  (5- 12d) 

n  =  1  ,2,3 
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First-Order  Solutions 

Two  of  the  first-order  constants  of  integration  are  evaluated  by 
applying  (5-2)  to  (4.6). 

7,(0)  =  xa  [ca  cos(ut)  -  Aa  sin(u)t)]  =  0,       (5.13) 

where 

*a  =  6ea/Bqla-  (5.1*0 

Since  (5.'3)  must  hold  for  all  values  of  time, 

Aa  =  Ca  =  0.  (5.15) 

The  other  two  constants  are   determined  by  applying  (5- 12a)  and  (5.15) 
to  (4.7). 

v,a(0)  =  Ba  cos(u-t)  +  Da  sin(ut)  =  (a/2)  sin(cot)  (5-16) 

Hence, 

Ba  =  0  (5.17a) 

Da  =  a/2.  (5.17b) 

The  first-order  variations  are  obtained  by  substituting  (5- '5)  and  (5-17) 
into  (4.6)  and  (4.7). 

T]a  =  xa  (a/2)  sin(z)  cos(cat  -  xgz)  (5-l8a) 
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where 


v,a  =  (a/2)  cos(z)  sin(ut  -  xai) ,  (5.18b) 


i=Bq]az.  (5.19) 


These  solutions  are  identical  with  the  small-signal  solutions  obtained 
by  Hahn  [l ]  and  Ramo  [2]. 

Second-Order  Solutions 

Before  obtaining  the  second-order  solutions,  it  is  necessary  to 
evaluate  Kj  through  Kg  as  given  in  Appendix  B.   The  values  in  Region 
A  are  obtained  by  substituting  (5-15)  and  (5-17)  Into  the  general 
express  ions . 

Rla  "  R3a  =  'Kka  =   R8a  =  °  <5-20a) 

K2a  =  l{cL/h)2  (5.20b) 

K   -  -  3  (aA)2  (5.20c) 

ba 

K6a  =  xa(aA)2  (5.20d) 

K7a  =  3xa(aA)2  (5.20e) 

Two  of  the  second-order  constants  of  integration  are  obtained  by  substi- 
tuting (5.2)  and  (5.20)  into  (4.18). 


Hence, 


54 


4a(0)  =  2  xa  [Ha  +   Va  ^/k)      0  +  3  9a)J  cos(2ut) 

-  2  xaJa  sln(2ut)  =  0.  (5.21) 

Ja  =  0  (5.22a) 

"a  =  "  Va^  °  '  Sa2)  (a/2)2 '  <5-22b) 


The  remaining  second-order  constants  are  determined  by  substituting 
(5.12b)  and  (5.20)  into  (4.19)  and  comparing  functions  of  time. 


*2a(0)  =  Ea  +  {<Z/h)2   '  (,/2)  {a/k)1  +   (VSa}  sin(2ut) 


Hence, 


+  [Ta/sa  -  («/4)2  -  (1  +  2  s  2)  ga  (a/4)2] 


cos(2ut)  =  -  (a/4)2  [1  -  cos(2wt)].       (5.23) 


Ea  =-(3/8)  (a/2)2  (5.24a) 


Ta  =  3  (a/4)2  saga  (1-2  sa2)  (5.24b) 


La  =  0.  (5.24c) 


In  the  discussion  of  boundary  conditions,  it  was  mentioned  that  the 
initial  polarization  (z  )  is  zero  provided  the  grids  are  closely  spaced. 
The  remaining  constant (Ha)  is  determined  by  applying  this  boundary 
condition, (5. 20),  (5-22),  and  (5.24)  to  (4.17). 
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Ha  =  0  (5.24d) 

The  second-order  current  is  found  by  substituting  (5.22)  and  (5.24)  Into 
(4.18). 

i2a  =  (1/2)  (a/2)2   ta2[l  +  3  ga  cos(2z)  -  4  ga  (1  -  sa2) 

•  cos(z/sa)]  cos  2(ut  -  xaz)  +  (3/2)  (a/2)2 
'    ^a  [sa  sfn(2i)  -  (1  "  2  sa2)  sin(z/sa)] 

•  sin  2(ut  -  xaz)  (5.25) 

The  second-order  velocity  is  found  by  substituting  (5.22)  and  (5.24) 
into  (4. 19). 

v2a  =  -  (a/4)2  [cos2  (z)  +  (1  -  3  ga  (1  -  2  sg2)  cos(z/sa) 
+  9*  0  +  2  s=.2)  cos(25))  cos  2(ut  -  X  5)  +  x  a 

a  a  a         a  a 

•((1  +  2  sa2)  sin(2l)  -  4  sa(l  -  sa2)  sin(z/sa)} 
•sin  2(ut  -  xaI)]  (5-26) 

Third-Order  Solutions 

In  order  to  evaluate  the  third-order  solutions,  the  parameters  Kg 
through  K24,  as  defined  in  Appendix  B,  must  be  evaluated  for  Region  A. 
Equations  (5-15),  (5-17),  (5-20),  (5.22),  and  (5-24)  are  used  to  evaluate 
these  parameters . 
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Kga  =   (9/4)    («/2)3   tasag. 


(5.27a) 


Kin     =    (3/4)    (a/2)3   x  g      [2   s      (b        -   s   l)    -    (    6   s      +    l)/(4  s   )]    (5.27b) 
10a  a   a  a        2a  a  a  a 


KHa  =    -    (3A)    {a/l)'   xaga    [2   sa    (bla    -   s/)    +   (6   sa    -    l)/(4   sa)] 


R]2a=    (15A)    (0=72)3    Satag£ 


K13a  =   Kl4a  =   Kl 5a  _   K16a  "  ° 


K17a  =   Kl8a  =   K19a  =   K20a  ~  ° 


2    .    „        2 


(5- 

27c) 

(5. 

27d) 

(5. 

27e) 

(5. 

270 

2    .    2 


K,.     =  (a/ky  ga   [54  s/  +  2  x/  +  14  t/  -   16  s/   t/  -  9J/2         (5-27g) 

Z  I  d  a  a  a  a  a  a 

iL_       -  2   (aA)3   b)a   (ga/sa)    [9(1-2   sa2)/4  +  xatab,a   (b^  -  2s/)] 


x22a 


a   '  a  a  2a    Ma  a 

(5.27h) 


K23a  =  -  2  (aA)3  b2a   (ga/sa)    [9(1-2  sa2)/4  +  xatabla(b3a  ♦  2  sa2)] 


*24a  "   <*A)3   9a    [9°   ^   "  9  +   l4  ta'   "  2  Xa2]/2 


(5.270 
(5-27J) 


Two  of  the  third-order  constants  of  integration  are  evaluated  by 
applying  (5.2)  to  (4.27). 

T3(0)  =  xa  G3(0)  cos(cot)  -  xa  F3(0)  sln(ut)  -  0        (5-28) 

Since  this  condition  must  hold  for  any  value  of  time, 


G3(0)  =  0 


(5.29a) 
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F3(0)  =  0.  (5.29b) 

Applying  (5-29)  to  (4. 24)  and  (4.25)  gives 


P"a  +  R10  sa2/b3a  "  Rl  1  sa2/b4a  '   R12/8  =  °  (5.30a) 


*a  +  R18  Sa2/b3a  *  R19  S^a  "  R20/8  =  °'  <5'30b> 


The  third-order  constants  Pa  and  Ra  are  found  by  transposing  (5-30)  and 
applying  (5- 27f ) . 


?a=  R12/8+sa2  [Rll/b4a  "  R10/b3a]  ^'^ 


Ra  =  0  (5.32) 


The  remaining  third-order  constants  are  determined  by  applying  (5- 12c) 
to  (4.28). 


V3a(0)  =  [F3a  (0)  +  Faa(0)  +  Fca(0)]  cos(ut) 


+  [G3a'(0)  +  Fba(0)  +  Fda(0)]  sin(u)t) 
=  3  (aA)3  sin(cjt),  (5-33) 


where  F3  (0)  implies  differentiation  of  Fo(z)  with  respect  to  z  and 
evaluation  at  z  =  0.   Applying  (5.27)  to  the  derivatives  of  (4.24) 
and  (4.25)  gives 
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F3a'(0)  =Qa  (5-34a) 


G~3a  (°)  =  S"a  +  (sabla/b3a)  R22a  "  R21a/2  "  (sab2a/b4a) R23a 

-  (3/8)  l2ka.  (5.3^+b) 

The  values  of  Fa(z)  through  F£j(z)  for  Region  A,  and  evaluated  at  z  =  0,  are 

Faa(0)  =  0  (5.35a) 

Fba(0)  =  -  9  (a/4)3  (5.35b) 

Fca(0)  =  0  (5.35c) 

Fda(0)  =  -  3  (a/4)3.  (5.35d) 

Substituting   (5.3*0    and   (5-35)    into   (5-33)  gives 


Qa   cos(tt)    +    [Sa   -   R2)a/2   +   (sab]a/b3a)    R22a   -    (sab2a/b4a) 

•K.,      -   (3/8)    JL.      -    12   (a/4)3]   sin(ut)   =  3    (a/4)3   sfn(ut). 
23a  24a 

(5.36) 
The  third-order  constants  are  evaluated  by  comparing  functions  of  time. 


Qa  -  0  (5.37) 

§a  =  *21a/2  "  ^Sabla/b3a)  *22a  +  (3/8>  R24a  +  (Sab2a/b4a>  R23a 

+  15  («/4)3  (5.38) 
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Since  the  third-order  constants  have  been  established,  the  third- 
order  solutions  for  the  current  and  velocity  can  be  evaluated. 

!3a  =  V*3  [A6  sin^)  +  A7  sin(b)az/sa)  +  A8  sln(b2az7sa) 

+  A9  sin(3z)  +  A]0  2  cos  (I)]  cos(ut  -  x  i)  -  x  «3 

a      a 

•  [A,  cos(z)  +  A2  cos(b)a2/sa)  +  A3  cos(b2az/sa)  +  A^ 

•  cos(3z)  +  A5  z  sin(I)]  sin(ut  -  x&z)  (5.39) 

v3a  "  a      CAn  z  cos(i)  +  A,2  sin(z)  +  A,3  s  in(b,az"/sa) 

+  A]Z+  sin(b2az/sa)  +  a15  sin(3z)]  cos(u)t  -  xaz)  +  a? 
•CA16  z  sin(z)  +  A,7  cos(z)  ♦  A]g  cos(b|ai/sa)  ♦  Al9 
•cos(b2ai/sa)  +  a20  cos(3i)J  sln(ut  -  xgi)  (5.40) 

The  constants  A,  through  A2Q  are  evaluated  in  Appendix  B. 


Results 

The  fundamental  component  of  current  is  found  by  adding  (5. 18)  and 
(5.39). 

!fa  "  ha  +  '3a  (5.M) 
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In  Fig.  5-'  the  amplitude  of  the  fundamental  frequency  current  is  plotted 
against  distance  along  the  beam  axis  (z)  for  three  values  of  a.      The 
small-signal  solution  (i.  )  is  also  shown  for  the  purpose  of  demon- 
strating saturation  effects.   Several  interesting  points  are  noted 
from  examination  of  this  figure.   These  are   separately  discussed  below. 
The  amplitude  of  the  fundamental  frequency  velocity,  along  with  the 
corresponding  small-signal  solution,  is  shown  in  Fig.  5-2. 

Saturation  of  fundamental  current 

Small-signal  theory  predicts  an  increase  in  the  amplitude  of  the 
fundamental  current,  which  is  linearly  proportional  to  d\    however,  it 
is  apparent  from  Fig.  5-1  that  the  amount  of  increase  becomes  pro- 
portionally smaller  as  a  increases.   This  saturation  effect  is  a  well- 
known  experimental  fact.   There  is  some  reason  to  doubt  that  the 
apparent  growth  in  current  with  distance  for  o.   =  0.3,  as  ascertained 

from  this  theory,  is  physically  real,  inasmuch  as  this  data  comes  from 

2 
a  region  where  overtaking  has  occurred.   Fig.  5-3  shows  the  effect  of 

saturation  as  maximum  current  is  plotted  against  ct.       It  should  be 

noted  that  the  apparent  Increase  in  the  slope  past  ct   =  0.3  is  well 

outside  the  range  where  this  third-order  approximation  is  valid. 

Locations  of  velocity  and  current  maxima  and  minima 

It  is  also  apparent  from  Fig.  5-1  and  Fig.  5.2  that  the  locations 
of  the  maximum  and  minimum  current  and  velocity  shift  as  &    is  increased. 
The  locations  of  the  maximum  and  minimum  values  of  current  and  velocity 


^The  concept  of  overtaking  is  discussed  in  Appendix  C. 
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Fig-  5.1  -  Fundamental  frequency  current  vs.  posit 
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Fig.  5.2  -  Fundamental  frequency  velocity  vs.  position. 
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Fig-  5.3  -  Maximum  fundamental  frequency  current  as  a  function  of  a. 
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as  a  function  of  a.  are  shown  in  Fig.  5-4.  These  shifts  from  the  iocations 
predicted  by  small-signal  theory  are  appreciable  even  for  a  moderate  value 
of  a.  such  as  0.2.  This  shift  in  position  has  been  experimentally  observed 
by  Mihran  [13]. 

Phase  shift 

Since  this  analysis  is  a  nonlinear  approach,  phase  shift  should  be 
expected.   The  fundamental  current  is  given  by  the  sum  of  (5-18)  and  ( 5 • 39) » 
and  may  be  expressed  in  the  form 


ir  =  lr   cosfcot  -  x   z  -  &,) . 
fa    fa  a      f 


(5-42) 


Examination  of  (5-39)  shows  the  presence  of  both  a  cos(u)t)  and  a  sin(u)t) 
component.  The  first-order  term,  as  given  by  (5-18),  has  only  a  cos(u)t) 
component.   Hence,  the  magnitude  of  the  fundamental  current  is  given  by 


',    =^(T,  +  T,  )2  +  (7.  )2 

fa      la    3ac       3as 


(5.^3) 


where  the  subscripts  £  and  s_  indicate  cos(o;t)  and  sln(cot)  coefficients, 
respectively.   The  phase  shift  is  given  by 


'fa 


=  tan' 


'3as 


la    3ac 


(5.44) 


The  phase  shift  increases  markedly  with  increasing  a,    as  indicated  in 
Fig.  5-5.   This  phase  shift  is  evaluated  at  the  point  of  maximum  current 


65 


100 


Fig.  5.4  -  Location  of  current  and  velocity  maxima  and 
minima  as  a  function  of  a. 
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Fig.  5.5  -  Phase  shift  at  the  point  of  maximum  current  as  a 
function  of  a. 
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for  each  value  of  a.      These  results  indicate  that  the  technique  of  ignor- 
ing the  out-of-phase  component,  as  employed  by  Paschke  [23,  2*+],  leads 
to  considerable  error,  particularly  for  larger  values  of  a. 

Second  harmonic  generation 

The  magnitude  of  second  harmonic  current  is  shown  in  Fig.  5.6. 
Examination  of  (5.25)  indicates  that  the  magnitude  of  second  harmonic 
current  increases  as  the  square  of  the  value  of  a.      Comparison  of  Fig.  5-1 
and  Fig.  5.6  would  tend  to  indicate  that  better  performance  could  be 
obtained  by  tuning  the  klystron  output  cavity  to  the  second  harmonic 
rather  than  the  fundamental.   This  conclusion  is  unwarranted,  however, 
since  the  fundamental  current  solution  is  based  on  terms  from  the  first- 
order  and  third-order  solutions,  with  the  third-order  term  contributing 
to  saturation.   On  the  other  hand,  the  second  harmonic  solution  shown 
neglects  an  additional  second  harmonic  term  which  would  be  contributed 
by  a  fourth-order  solution,  and  therefore  fails  to  show  saturation  effects 
The  second  harmonic  solution  obtained  by  this  third-order  approximation 
is  analogous  to  the  fundamental  frequency  solution  obtained  by  small - 
signal  analysis. 


68 


K  =  1.0 

2/3 


2.0 


1.6  -" 


1.2 


0.8 


0.4  ■- 


Fig.  5.6  -  Second  harmonic  current  vs.  position. 


CHAPTER  6 

PARTICULAR  SOLUTIONS  FOR  REGION  OF  SECOND  DRIFT  TUBE 

Synops  is 

In  this  chapter  methods  are  investigated  for  amplifying  the  current 
beyond  that  obtained  in  the  simple  two-cavity  klystron  discussed  in  the 
previous  chapter.   The  model  for  the  space-charge-wave  amplifier  is 
identical  to  a  two-cavity  klystron,  except  the  drift  tube  is  divided 
into  sections,   each  of  which  may  have  a  different  potential  and  a 
different  diameter.   The  additional  amplification  results  from  changes 
in  the  drift  tube  potential  and  radius.   Both  of  these  changes  effect 
the  value  of  the  plasma  reduction  factor  (R]).   This  amplifier  represents 
a  combination  of  the  velocity-jump  amplifier  proposed  by  Tien  and  Field 
[3],  and  the  rippled-wall  amplifier  investigated  by  Blrdsall  and  Whinnery 
[k].      A  third  method  for  obtaining  amplification  is  to  vary  the  beam 
diameter;  however,  this  technique  is  not  considered  because  of  the  addi- 
tional complexity  introduced  both  in  the  theoretical  analysis  and  in 
construction  of  a  working  unit.   The  analysis  includes  the  effect  of 
placing  the  transition   in  either  the  vicinity  of  maximum  velocity  modu- 
lation, or  In  the  vicinity  of  minimum  velocity  modulation.   The  investi- 
gation also  includes  a  search  for  the  optimum  location  of  the  transition, 
in  so  far  as  concerns  the  achievement  of  maximum  current  modulation  at 
some  point  further  down  the  beam.   The  solutions  for  current  and  velocity 
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are  obtained  from  the  general  solutions  in  Chapter  k,    by  applying  boundary 
conditions  at  the  transition.   All  data  for  the  figures  in  this  cl apter 
are  obtained  from  computations  performed  by  the  IBM  709  computer.   The 
data  for  the  input  to  the  transition  are  obtained  from  the  solutions  for 
Region  A,  using  the  beam  parameters  specified  in  Chapter  S- 


Boundary  Conditions 

The  modulated  electron  beam  drifts  under  the  influence  of  space- 
charge  forces  until  it  reaches  the  position  z  =  Z  ,  at  which  point  the 
second  drift  tube  is  positioned.   The  drift-tube  potential  changes  from 

V  to  V,  at  the  junction  of  the  two  tubes.   From  the  conservation  of 
a     b 

energy, 

(m/2)  (v,  2  -  v  2)  =  lei  (V.  -  V  )  , 
v  b     a  '       v  b    a' 


V  =  2  ^  <Vb  "  Va>  +  va2-  (6-'> 


Using  the  relation 


v02  =  2  n  V,  (6.2) 


■*To  help  insure  that  the  computer  program  for  Region  B  was  correct, 
computations  were  performed  without  making  any  change  in  either  drift- 
tube  potential  or  plasma  reduction  factor.   The  results  of  these  compu- 
tations were  compared  with  the  solutions  for  Region  A  extended  through 
the  distance  covered  by  Region  B.   The  results  were  identical  to  three 
significant  figures.   The  solutions  obtained  for  Region  A  were  compared 
with  the  solutions  obtained  by  Engler  [26]. 
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(6.1)  may  be  expressed  as 


V,  =^V  ,  2  -  V   2  +  V  2  ,  (6.3) 

h       oh       na       a  v   -^ ' 


b     ob     oa     a 


v.  =  v  .  ^1  +  (v  2  -  v   2)/v   2  .  (6.4) 

b    ob      v  a     oa  '      ob  v   ' 

Expanding  (6.4)  in  a  binomial  series  gives 

vb  =  vob  C1  +  y/2  -  y2/8  +  y3/'6],  (6.5a) 

where 

y  =  (va2  -  voa2)/vob2-  (6-5b) 

The  velocity  (va)  may  be  written  as  a  power  series  in  a   according  to 
the  method  of  successive  approximations. 

va  =  voa  +  vla  +  v2a  +  v3a>  (6-6) 

where  the  various  order  terms  (v  )  are  proportional  to  an .      Substitution 
of  (6.6)  into  (6.5b)  yields 

y  =  £2  vna  (vl,  +  v9.  +  viJ  +  vi.2  +  v9=2  +  v-    J 


oa  v  la    2a    3a     la     2a     3a 

',   v0  +  v.   v,   +  v0   v,  )]/v  ,i 
la  2a    la  3a    2a  3a'   ob 


+  2  (vl.  v?,  +  v..  v^  +  v?a  v^)]/vJ'  (6-7) 


It  is  assumed  that  the  a-c  terms  are  all  considerably  smaller  than  the 
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d-c  terms.    Consequently,  (6.7)  may  be  approximated  by 


V  -  2  voa  <vla  +  v2a  +  v3a)/vob*-  (6-8) 


Equation  (6.8)  is  substituted  into  (6.5a),  and  the  approximation  that 
vna/voa  «    1  is  applied.   The  resulting  velocity  in  Region  B  is 

vb  =  vob  []    +   voa  (vla  +  v2a  +  ^'^ob*2!-      (6"9) 

The  method  of  successive  approximations  is  applied  to  v.  .   Comparison 
of  terms  in  (6.9)  yields 


v   (Z  )  =n/v  /V.  v   (Z  ).  (6.10) 

nbx  a'     a   b  nav  a'      ,   <•>   -, 

n  =  1 ,  2,  3 


Essentially  this  is  the  small-signal  assumption  used  in  the  analysis  of 
the  velocity-jump  amplifier  by  Tien  and  Field;  however,  the  solution  in 
Region  B  contains  second-order  and  third-order  terms  In  addition  to  the 
small-signal  solution. 

The  second  boundary  condition  is  the  continuity  of  current.   Hutter 
[ll]   shows  that  the  current  is  continuous  across  a  velocity  jump  by 
using  kinematic  considerations,  and  alternately  as  a  limiting  case  of 
the  Llewellyn-Peterson  equations. 


For  a  =   0.2  the  largest  a-c  velocity  term  is  approximately  10  per 
cent  of  the  d-c  velocity. 
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Consequently, 

i    .     (Z    )    =    i        (Z   ).  (6.11) 

nb    v    a'  na    v    a7  v  ' 

n  =    1 ,    2,    3 

First-Order   Solutions 
The   first-order   current    in   Region   B    is    found   from   (4.6). 

Tj.    =   BEb    [Cb  cos(BQ.lb   i)    +  Db   sin(BQlb   z)  ]    cos(cot    -   BEb   z) 

-   BEb    [A\    cos(BQlb   z)    +  b\    sin(BQJb   z)  ]    sin(ut    -   BEb   z) 

(6.12) 
The   first-order   current    in   Region  A    is   obtained   from   (5>l8a),    and    is 
evaluated  at    the  boundary   z  =   Za- 


i]a(Za)    =    (a/2)    xa   sin(Za)    cos(oot    -  xq   Za)  (6.13) 


T\-jo   equations  for  the  first-order  constants  are  determined  by  applying 
the  continuity  of  current  condition  to  (6.12)  and  (6.13).   The  two 
equations  are  obtained  by  comparing  coefficients  of  cos(u)t)  and  sin(o)t) 


(a/2)    xa  sin(Za)  cos(xa  Za)  = 


SE.  [C,  cos(BQl.  Z  )  +  6.  sin(BQ.l.  Z  )]  cos(BE.  Z) 

b   b        b  a     b        oa        Da 


+  BEL    [AL    cos(BQl      Z   )    +   BL    sin(BQl      Z   )]    sin(BEL   Z   ) 
bb  bab  ba  ba 


(6.14a) 


lh 


(a/2)    xa   sin(Za)    sin(xa   Za)   = 

BEb    [Cb   cos(BQlb   Za)    +   Db   sin(BOJb   Za) ]    sin(BEb   Za) 

-   BEL    [A.    cos(BQl.    Z   )    +   B,    s  Jn(BQ.l.    Z    )]    cos(BE     Z   ) 

bb  b      a  d  Da  Da 

(6.14b) 

The  other  two  equations  needed  to  evaluate  the  four  constants  are 
derived  from  the  velocity  boundary  condition.   The  first-order  velocity 
in  Region  A  is  obtained  from  (5.1 8b) . 


v   (Z  )  =  (a/2)  cos(Z  )  sin(ut  -  x  Z  )        (6.15) 
la  a  a  a  a 


The  first-order  velocity  in  Region  B  is  found  from  (4.7) 


v,.(z)    =n/v.  /V      BOJ.     [B.    cos(BQJ,    z)    -  A,    sin(BQJ      z)  ] 
lb  babb  b  d  d 

•  cos(ut  -  BEb  z)  +  J\/V~a   BQ1  b  [5b  cos(BQlb  z) 


CL  sin(BQl   z)]  sin(ut  -  BE.  z)  (6.16) 

b        b  b 


The  first-order  velocity  (v)b)  is  normalized  with  respect  to  the  d-c 

velocity  (v   )  in  Region  A,  rather  than  with  respect  to  the  d-c 
'    oa 

velocity  (v  b)  in  Region  B.   This  is  done  in  order  that  the  two  normalized 
first-order  velocity  terms,  vQa  and  v  .,  may  be  compared  directly.   The 
boundary  condition  (6. 10)  is  applied  to  (6.15)  and  (6.16).   The  remaining 
two  equations  needed  to  evaluate  the  first-order  constants  are  obtained 
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by  comparing  coefficients  of  cos(wt)  and  sin(cot). 

-  (a/2)  (Va/Vb)  cos(Za)  sin(xa  Za)  = 

BQlb    [Bb  cos(BQlb  Za)    -  Ab   sin(B0Jb  Za) ]    cos(BEb  Za) 

-   B0Jb    [Db   cos(BQlb   Za)    -   Cb   sin(BQlb   Za)  ]    sin(BEb  Zg) 

(6.17a) 

(a/2)    (Va/Vb)    cos(Za)    cos(xa   Za)   = 

BQlb    [Bb   cos(BQlb   Za)    -  Ab   sin(BQlb  Zg) ]    sin(BEb  Za) 


+   BQ.1.     [6,     cos(B0J,     Z   )    -   C,     sin(BQI,    Z   )]    cos(BE,    Z   ) 
bo  bab  ba  ba 


(6.17b) 


The  first-order  solutions  are  established  once  the  constants  A,  ,  B.  ,  C,  , 

b   b   b 

and  Db  are  evaluated  and  substituted  into  (6.10)  and  (6.16).   Cramer's 

5 
Rule   is  used  to  evaluate  the  constants. 


Second-Order  Solutions 
The  second-order  current  in  Region  B  is  obtained  from  (^4.18) 

i2b(z)  =  2  BEb  tY]b(z)  +  Mb  cos(BQ2b  z)  +  Lb  sin(BQ2b  z) ] 

•  cos  2(«ut  -  BEb  z)-  2  BEb[Y2b(z)  +  Jb  cos(BQ2b  z) 


^See  Appendix  D  for  a  discussion  of  the  solutions  for  these  constants 
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+  ?b   sin(BQ2b   z)  ]    sin  2(ut    -   BE      z) ,  (6.18) 

where 

Ylb(*>  =  Sb  [R6b  +  9b  S  Cos(2  BQ,b  'z)    +  %   *8b  s'n(2BQ'b  z)]/Ba2b 

(6.19) 

Y2b(i)   =   Sb    [R3b  +  9b   \b  C°S(2   BQ'b   Z_)    +  9b   R5b   Sin(2   BQ'b   ;)^Bd2bi 

(6.20) 

The  second-order  current  in  Region  A  is  obtained  from  (5.25),  and  evalu- 
ated at  the  transition  (Za). 

»ta<ZJ  "  Yla<Za>  cos  *(«*  "  xa  Za>  +  Y2a<Za) 

•  sin  2(ut  -  xa  Za)  ,  (6.21) 

where 

Y.  (Z  )  =  (1/2)  (a/2)2  t  2  [1  +  3  g  cos(2  Z  ) 

13    3  3  3  3 

-  k   ga  (I  -  sa2)  cos(Za/sa)]  (6.22) 

W  =  (3/2)  (a/2)2  'a  9a  [sa  Sin(2  Za> 

-  (1  -  2  sa2)  sin(Za/sa)].  (6.23) 

Equation  (6.18)  is  equated  to  (6.21),  using  the  boundary  condition  for 
continuity  of  current  given  by  (6.11),  and  coefficients  of  cos(cot)  and 
sin(d;t)  are  compared. 
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Yla(Za)    cos(2   xa   Za)    -   Y2a(Za)    sin(2   xa    Za)    = 

2BE      [Y,L(Z    )    +   K     cos(BQ2      Z   )    +   L,     sin(BQ2      Z    )] 
b        1  b     a  b  ba  b  t>     a 

cos(2    BEk  Zj    +  2    BEk    [Y,k(Zj    +   Jk  cos(BQ2k  Zfl) 

+   T      sin(BQ2      Z   )]    sin(2BE     Z   )  (6.24a) 

b  b     a  b     a 


Y,a(Za>    sin<2  Xa   Za>    +  Y2a<Za>    cos<2  Xa   Za>    = 

2    BEb    [Ylb(Za)    +  Sb  cos(B(>2b  Za>    +   [b   sin(B(>2b   Za)] 

•  sin(2  BEk  Za)  -  2  BEk  [Y2b(Za)  +  Jk  co5(BQ2k  Za) 

+  T   sin(BQ2k  Za)J  cos(2  BEb  Zfl)  (6.24b) 

The  second-order  velocity  in  Region  B  is  obtained  from  (4.19),  and 
normalized  with  respect  to  voa. 


v2b<z>  ='/VVa  [Eb  +  BEb  B^b  P10b/2  +  (B(i,b)   P2b/Z+ 

-  (K2b/6)  cos(2  BQlb  z)  +  (K]b/6)  sin(2  BQ1 b  z) ] 


Wv./V   [BQ2b  (T  cos(B0.2k  5)  -  J"b  sin(B0_2b  z) } 


+  Y   (z)]  cos  2(ut  -  BEb  z)  +  >/  Vb/Vg  [BQ2b  {Lb 

•  cos  (BQ_2b  z)  -  Mb  sin(B0.2b  z) }  +  Y^z)] 

•  sin  2(ut  -  BEb  z),  (6.25) 


78 


where 

y,l>)  =  "  R^/xu  +  (sL/3)  ('  +  2  SL2)   "^  ,  cos<2  *Q-K   ^ 

3b        6b  b     b  b     5b         b 

-  Kkh   sln(2  BQlb  z)  3  (6.26) 

Y4b(z)  =  K3b/xb  +  (gb/3)  d+2  sb2)  [Rgb  cos(2  BQl b  z) 

-  K?b  sln(2  BQlb  z) ]  .  (6.27) 


The  second-order  velocity  in  Region  A  is  found  from  (5-26),  and  evaluated 
at  the  point  of  discontinuity  (Z  ). 

*?JZJ  ■  -  <aA>2  [cos2  (z.)  +  WZJ  ccs  2(u:t  "  x*  ZJ 

+  Y4a(Za)    sin  2(ut    -  xa   Za)],  (6.28) 

where 


Y*a<zJ    =    1    +  9,    (1    +  2   s   Z)    cos(2   Z  )    -  3  g      (1    -  2   s/)    cos(Z   /s    ) 
ja     a  a  a  aa  d  a     a 

(6.29a) 

WZJ   =  x*  ^      (1    +  2   s   2)    sin(2  Z   )    -  4  g     t      (1    -  s   2)    sin(Z  /s   ). 
4a     a  a     a  a  a  aa  a  aa 

(6.29b) 


Equations  (6.25)  and  (6.28)  are  compared  at  z  =  Za  by  using  the  boundary 
condition  of  (6.10).   Three  equations  for  the  second  order  constants  are 
obtained  by  comparing  the  coefficients  of  cos(ut),  and  sin(u)t),  and  the 
time  independent  terms. 

-  {Ci/k)1    (Va/Vb)  cos2  (Za)  =  Eb  +  (1/2)  BEb  BQJ  b  P|ob  +  (1/4) 
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(BQlb)2   P2b   -    (1/6)    K2b  cos(2    BQlb  Zg)    +   (1/6)    K,b 

sin(2   BQlb   Za)  (6.30) 


-(a/k)2    (Va/Vb)    [Y3a(Za)    cos(2  xg   Zfl)    -   \atta)    sin(2  xa   Zfl)  ] 
=    [BQ2b    (?b  cos(BC)2b  Zg)    -   Jfa   sin(BQ2b  Za) }    +  Y3b(Za)] 

•  cos(2    BEb  Za)    -    [BQ2b    {    Lb   cos(BQ2b  Za)    -  Mb  sin(B0_2b  Za) } 
+  Y4b(Za)  ]  sin(2   BEb  Za)  (6.31a) 

-(«A)2    (Va/V„)     [Y3a(Za)    sin(2  xa   Za)    +  Y4a(Za)    cos(2   xa   Za)  ] 
=    [BQ2b    [7b  cos(BQ2b  Za)    -   \   sin(BQ2b  Zg)}    ♦  Y3b(Za)] 

•  sin(2   BEb  Zg)    +    [B0_2b    (Lb  cos(B0_2b   Za)    -   Mfa  sin(B0_2b  Za) ) 
+  Y4b(Za)]    cos(2   BEb  Za)  (6.31b) 

The   second-order   constant    Eb  may   be  obtained   from   (6.30). 


Eb  «=    -    (OA)2    (Va/Vb)    cos2    (Za)    -    (1/2)    BEb   BQlb   P,ob    -    (I  A) 
"    (B(ilb)2   P2b  +   C/6)    K2b  cos(2   BQlb   Za)    -    (1/6)    K]b 

•    sin(2    BQJ.    Z    )  (6.32) 

b     a 


Equations  (6.24)  and  (6.31)  represent  a  set  of  four  equations  in  the 
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four   second-order   constants.      The   solutions    to  this   set  of  equations   are 
also  discussed    in  Appendix   D. 


Third-Order   Solutions 


The   third-order   current    in   Region   B    is   obtained   from   (^.27) 


i      (z)    =   BE      [R     cos(B01      z)    +  s"      sin(B01      z)    +  Yr.(z)] 
3b  b       b  b  b  b  5b 

•    cos(oJt    -   BEb   z)    -   BEb    [Pb  cos(BQlb   z)    +  Qb 


sin(BQl      i)    +Y6b(i)J   sln(wt   -   BEb   z) ,  (6.33) 


where 


Y      (z)   =    (z/2)    [R  sin(B0.1b   z)    -   K2)b   cos(B0Jb   z)  ]   +   (s^b^) 

•  [K,8b  cos(BQ2b   b]b   z)    +  R22b   sin(BQ2b  b]b   z) ]/BQ2b 

"    (sb/bW    [i<19b  cos(B02b  b2b   z)    +   K23b   sin(BQ2b   b2b   z)  ] 

/BQ2b   -    [K20b  cos(3    BQlb   z)    +  K2kh   sin(3    B0Jb  z)]/(8  B0Jb) 

(6.34) 


Y6b(z)    -    (z/2)     [R^  sin(BQlb   I)    -   K^  cos(B0Jb   z)  ]   +   (s^) 
•  [K)0b  cos(BQ2b   b|b  I)    *  K]kb   sin(BQ2b   b]b   z) ]/B0_2b 
-    (sb/b4b)    [R)lb  cos(BQ2b   b2b   i)    +  K)5b   sin(B0_2b  b2b   i) ]/BQ2 
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-  [K|2b  cos(3  BQlb  z)  +  K|6b  sin(3  BQ1   z) ]/(8  BQ1  )  . 

(6.35) 

The  third-order  current  in  Region  A  is  obtained  from  (5-39)  and  evaluated 
at  the  transition  (Za) . 


,3a(Za)  =  Y5a(Za}  COs(wt  "  Xa  Za>  "  Y6a(Za}  sin(wt  "  *a  Za) • 

(6.36) 

where 

Y5a(Za)  =  xa  ^    [A6  sin(Za>  +  A7  sin<bla  Za/sa)  +  A8  sin(b2a  Za/sa) 

+  A9  sin(3  Za)  +  A,0  Za  cos(Za)]  (6.37) 

Y6a(Za)  "  xa  ^    [a1  cos(Za)  +  A2  cos(bla  Za/sa>  +  A3  Cos(b2a  Za/sa> 
+  A^  cos(3  Za)  +  A5  Za  sin(Za)].  (6.38) 


Equations  (6.33)  and  (6.36)  are  compared  at  z  =  Za  by  using  the  continuity 
of  current  boundary  condition.   Equating  coefficients  of  cos(u)t)  and  sin(ut) 
yields  two  of  the  equations  needed  to  evaluate  the  third-order  constants. 


Yr  (Z  )  cos(x  Z  )  +  Y.  (Z  )  sin(x  Z  )  = 
5a  a      a  a     6a  a      a  a 


BE.  [R.  cos(BQl.  1)    +   S.  sin(BQj.  Z)    +   Yc.  (Z J] 


cos(BEb  Za)  +  BEb  [Pb  cos(BQlb  Za)  +  Qb 

sin(BQlb  Za)  +  Y6b(Za)]  sin(BEb  Za)  (6.39a) 
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VZa)    sin(xa   Za>    "   Y6a<Za)    cos(xa   Za>   = 

BEb    [Rb   cos(BQlb   Za)    +   §b   sin(BQlb   Za)    +  Y5b(Za)  ] 

•    sin(BE.    Z   )    -   BE,     [P,    cos(BO_l,    Z   )    +  Q_,     sin(BQ1L   Z   ) 
b      a  bb  ba  b  b     a 

+  Y6b(Za)]    cos(BEb   Za)  (6.39b) 

The   third-order   velocity    in   Region  A   is   obtained   from   (5.40)    and 
evaluated  at    z  =   Za. 


Sa(Za)    =  W    COS(uJt    "  *a   Za}    +  Y8a(Za)    Sin(U*    "  *a   ZaM6^0> 


where 


Y7a(Za)  =  ^  [A11  Za  COS(Za}  +  A12  »",<V  +  A.3  Sin(bla  Za/sa) 

+  A..  sin(b„  Z  /s  )  +  Air  sin(3  Z  )]  (6.4la) 

IH      L°      a  a     15       3 


YQ  (Z  )  =  a3  [A,.  Z   sin(Z  )  +  A.,  cos(Z  )  +  A.Q  cos(b.   Z  /s  ) 
oa  a        16  a      a     17      a     18      la  a  a 

+  A   cos(b2a  Za/sa)  +  A2Q  cos(3  Za)J.  (6.41b) 


The  third-order  velocity  in  Region  B  is  obtained  from  (4.28).   The  third- 
order  velocity  is  normalized  with  respect  to  vQa  rather  than  v  ,.   This 


Is  accomplished  by  multiplying  (4.28)  by  \l  V. /V_ 


v3b(z)  =  ^Vb/Va  [BQlb  Qb  cos(B0Jb  z)  -  BQlfa  p"b  sin(BQJb  z) 
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+  Y      (2)    ]   cos(wt    -    BE      i)    +    4TaT    [BQ1      Sl    cos(BQ1      z) 
1  o  babbb 

-   BQlb    Rb   sin(BQlb    z)    +   Y^)]    sfn^t    -    BEb   z)  ,  (6.42) 

where 

Y?b(i)    =   BQ)b    (z/2)    [R9b   cos(BQlb   z)    +  R^   sin(BOJb   z)  ] 

+  (b2b  W  [Rllb  S''n^\  b2b  i}  "  Rl5b  COS^2b  b2b  Z")] 
+  <blb  sb/b3b)  [Rl4b  cos(BQ2b  b,b  z)  -  K]ob  sin(BQ2b  b,b  \) ] 
+   (3/8)    [K,2b   s!n(3    BQlb   i)    -   R|6b   cos(3    BQlb   z)  ]   +   F ^{z) 

+  Fcb^)    +   CKgb/2)    sin(BQlb   z)    -    (K,3b/2)    cos(BQJb   z) 

(6.43a) 

Y8b(z")   =   BQlb    (z/2)    [R2]b   sfn(BQlb   z)    +   K??b   cos(BQlb   i)  ] 

+  (blb  Vb3b)  [R22b  cos(BQ2b  b,b  l)  -  R,8b  sin(BQ2b  blb  i)] 
+  <b2b  sb/b4b)  [Rl9b  sIn(BQ2b  b2b  z)  -  R23b  cos(BQ2b  b2b  z)  ] 
+   (3/8)     [R20b   sin(3    BQlb   z)    -   R2Afb  cos(3    BQJb   z)  ]   +  Fbb(z) 

+  Fdb(i)    +   (R!7b/2)    sfn(BQIb   Z)    "    (K2,b/2)    cos(BQlb   z") . 

(6.43b) 

The  remaining  pair  of  equations  required  to  evaluate  the  third-order 
constants  are  obtained  by  comparing  (6.40)  and  (6.42)  at  5  =  Za,  using 
the  boundary  condition  given  by  (6.10). 
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<Vvb>  [Y7a  <Za)  cos(*a  Za>  "  Y8a  (Za>  sin<xa  Za>  3  = 
BQlb  [$b  cos(BQlb  Za)  -  Pb  sin(BQlb  Za) ]  cos(BEb  Za) 

-   BQ1,     [S.    cos(BQl      Z   )    -   RL    sin(BQl      Z   )]    sin(BE      Z    ) 
b        b  bab  ba  ba 

♦  Y?b  (Za)  cos(BEb  Za)  -  Ygb  (Za)  sin(BEb  Za)  (6.44a) 

(VVb)  £Y7a  (Za)  sin<xa  Za)  +  Y8a  (Za>  cos(xa  Za>  ^  = 

BQ1.  [Q.  cos(BQl.  Z  )  -  P.  sin(BOJ,  Z  )]  sin(BE,  Z  ) 

u   u        d  a     d        Da         ba 

+  BQlb  [§b  cos(BQlb  Z  )  -  Rb  sin(BQlb  Z  )]  cos(BEb  Z  ) 

+  Y?b  (Za)  sin(BEb  Za)  +  Ygb  (Zg)  cos(BEb  Za)  (6.44b) 

Equations  (6.39)  and  (6.44)  comprise  a  set  of  four  equations  in  the  four 
third-order  constants.   The  solutions  for  this  set  of  equations  are  dis- 
cussed in  Appendix  D.   The  third-order  solutions,  given  by  (6.33)  and 
(6.42),  are  complete  once  these  four  constants  have  been  evaluated. 


Transition  Near  a  Velocity  Maximum 

According  to  the  small-signal  results  of  (2.16),  the  best  amplifi- 
cation for  a  transition  near  a  velocity  maximum  is  achieved  if  the  drift- 
tube  voltage  (V.)  is  reduced,  and  the  plasma  reduction  factor  (R,, )  is 
also  reduced  by  decreasing  the  drift-tube  radius  (a)  in  Region  B.   These 
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parameters  are  changed  In  relation  to  the  values  they  possessed  in  Region 
A.   The  conclusions  drawn  from  the  small-signal  analysis  are  that  best 
amplification  is  obtained  if  the  voltage  and  radius  are  changed  at  the 
exact  location  of  maximum  velocity  variation,  and  that  the  amplification 
obtained  Is  independent  of  a,    i.e.,  the  amplifier  does  not  saturate. 

The  effect  of  varying  the  drift-tube  radius  is  demonstrated  in 
Fig.  6.1  for  a  fixed  voltage  transition,  V.  =  Va/2.   It  is  seen  from 
Fig.  6.1,  that  a  10  per  cent  Increase  in  amplification  is  achieved  if 
the  b/a  ratio  is  changed  from  0.667  to  0.91  in  Region  B;  whereas,  the 
gain  drops  about  5  per  cent  if  the  ratio  is  changed  to  0.5  in  Region  B. 
Comparison  with  the  solution  for  Region  A,  which  corresponds  to  the 
conventional  two-cavity  klystron,  shows  that  a  25  per  cent  increase  in 
amplification  is  achieved  by  changing  the  b/a  value  to  0.9'  and  reducing 
the  potential  by  a  factor  of  two  in  Region  B.   The  major  contribution  to 
the  increased  amplification,  however,  comes  from  the  change  in  drift-tube 
voltage,  rather  than  the  change  in  drift-tube  radius. 

The  effect  of  varying  the  drift-tube  voltage  is  demonstrated  in 
Fig.  6.2  for  a  b/a  transition  from  0.667  to  0.91-   The  amplification  in- 
creases as  the  drift-tube  potential  in  Region  B  decreases.   The  curve  for 
V.  =  V  / 1 0  also  demonstrates  another  Important  consideration,  namely  the 
location  of  the  output  cavity.   If  the  position  of  the  second  cavity  is 
shifted  by  five  degrees  in  relation  to  its  optimum  value,  the  amplifi- 
cation Is  reduced  to  such  an  extent  that  a  much  smaller  voltage  transition 
would  yield  better  amplification.   The  actual  shape  of  the  curve  for  a  ten 
to  one  change  in  potential  Is  somewhat  In  doubt,  as  overtaking  occurs  over 
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Fig.  6.1  -  Current  in  Region  B  as  a  function  of  the  drift-tube  radius 
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Fig.  6.2  -  Current  in  Region  B  as  a  function  of  the  drift-tube  voltage. 
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a  considerable  portion  of  the  cycle.   It  is  also  evident  from  examination 
of  Fig.  6.2  that  the  optimum  length  of  the  second  drift  tube  is  a  definite 
function  of  the  voltage  ratio;  whereas,  the  curves  in  Fig.  6*1  show  that 
the  change  in  optimum  length  with  plasma  reduction  factor  is  negligible, 

Small-signal  theory  predicts  that  the  location  of  the  velocity-jump 
must  be  at  the  point  of  maximum  velocity  variation.   In  order  to  investi- 
gate how  the  predicted  amplification  is  affected  by  the  nonlinear  theory, 
the  location  of  the  transition  is  varied.   The  maximum  current  obtained 
at  some  point  in  Region  B  for  a  particular  transition  location  is  plotted 
in  Fig.  6.3  against  the  location  of  the  transition.   The  optimum  location 
for  a  =   0. 1  is  actually  about  twenty-two  degrees  closer  to  the  input 
cavity  than  the  point  where  maximum  velocity  variation  occurs.   The  op- 
timum amplification  is  about  5  per  cent  above  that  obtained  by  locating 
the  transition  at  the  point  of  maximum  velocity  variation.   It  should  be 
noted  that  the  position  of  the  cavity  is  not  too  critical,  in  that  a 
five  degree  shift  in  location  from  the  optimum  position  results  in  less 
than  a  1  per  cent  loss  in  amplification  for  this  example.   Thus  the 
contributions  to  the  over-all  amplification  due  to  changes  in  the  drift- 
tube  radius  or  location  of  the  transition  are  negligible,  at  least  for 
smal  1  values  of  cc. 

The  reason  for  this  increased  current  is  the  fact  that  while  the 
velocity  fluctuations  are  decreasing  as  the  transition  is  moved  toward 
the  input,  the  current  fluctuations  are    increasing.   Since  the  current 
is  continuous  at  the  transition,  the  current  excitation  in  Region  B 
builds  from  the  current  amplitude  at  the  transition.   Thus  a  compromise 
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Fig.  6.3  -  Maximum  current  in  Region  B  as  a  function  of  the  location 
of  the  first  transitu  n. 
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is  made  between  maximum  velocity  excitation  and  virtually  no  initial 
current,  and  the  case  where  little  velocity  excitation  is  achieved,  but 
a  large  current  is  present  at  the  transition.   In  the  small-signal  theory 
the  location  of  maximum  velocity  and  minimum  current  coincide,  but  this 
does  not  hold  for  the  nonlinear  analysis,  as  was  previously  shown  in 
Fig.  5.**. 

The  fact  that  the  velocity  at  the  transition  is  less  than  the  velo- 
city at  the  entrance  in  Region  A  indicates  that  this  amplifier  design 
might  not  be  correct.   The  curve  labeled  a     =  0.2  in  Fig.  6.1  supports 
this  conclusion.   This  curve  corresponds  to  the  current  obtainable  from 
a  conventional  two-cavity  klystron  with  the  same  signal  voltage  as  the 
space-charge-wave  amplifier,  but  employing  the  beam  characteristics  of 
Region  B,  namely  the  lower  d-c  potential  which  yields  a  beam  of  higher 
perveance.   To  exceed  the  current  output  of  this  klystron,  it  is  necessary 
to  use  an  amplifier  which  has  the  same  initial  beam  characteristics,  then 
employ  an  increase  in  drift-tube  potential  at  the  first  transition,  and 
finally  return  to  the  original  potential  at  the  second  transition.   The 
first  transition  is  made  near  a  velocity  minimum,  and  the  second  near  a 
velocity  maximum.   This  type  of  velocity-jump  amplifier  was  discussed  on 
a  linear  basis  in  Chapter  2,  and  the  small-signal  results  are  shown  in 
Fig.  2.3.   In  the  next  section  this  amplification  process  Is  investi- 
gated on  a  nonlinear  basis,  where  it  is  shown  that  amplification  is 
possible,  at  least  for  small  values  of  a. 
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Transition  Near  a  Velocity  Minimum 

The  smal 1 -signal  solutions  for  a  transition  at  a  velocity  minimum 
are  given  in  (2.9).   It  would  appear  that  the  variation  in  reduction 
factor  is  just  as  important  a  consideration  as  changing  the  potential 
because  of  the  four  to  one  ratio  in  the  powers  to  which  these  parameters 
are  varied.   Since  the  current  amplification  is  relatively  insensitive 
to  a  voltage  transition  at  a  velocity  minimum,  a  second  transition  is 
needed  in  the  vicinity  of  a  velocity  maximum.   As  shown  in  the  previous 
section,  the  best  amplification  is  achieved  at  the  velocity  maximum  if 
the  drift-tube  potential  is  reduced,  and  the  reduction  factor  is  decreas- 
ed by  reducing  the  drift-tube  radius.   Hence,  the  voltage  and  radius  of 
the  drift  tube  must  be  increased  at  the  first  transition,  then  returned 
to  their  original  values  at  the  second  transition. 

The  effect  of  varying  the  drift-tube  voltage  is  shown  in  Fig.  6.4. 
It  is  apparent  that  the  velocity  variations  are  relatively  insensitive 
to  the  change  in  potential;  however,  the  velocity  definitely  decreases 
as  the  drift-tube  voltage  in  Region  B  increases.   It  is  also  apparent 
that  the  optimum  location  of  the  second  transition  is  a  definite  function 
of  the  voltage  ratio.   The  fact  that  the  variations  are  relatively  In- 
sensitive to  changes  in  potential  follows  from  (2.9b).   The  ratio  of 
maximum  velocity  In  Region  B  to  maximum  velocity  in  Region  A  is  only  pro- 
portional to  the  one-fourth  power  of  the  voltage  ratio.   However,  it  also 
appears  at  first  glance  that  the  velocity  should  increase  with  increasing 
values  of  V^,  rather  than  decrease.   The  explanation  is  found  by  examining 
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the  plasma  reduction-factor  curves  shown  in  Fig.  2.1.   As  the  voltage 
increases  the  normalized  beam  radius  (6b)  decreases,  thus  decreasing 
the  value  of  the  reduction  factor.   Since  the  ratio  of  reduction  factors 
is  raised  to  the  first  power,  a  small  gain  due  to  the  increase  in  V^  is 
more  than  compensated  by  the  decrease  in  Ri.  . 

The  effect  of  varying  the  drift-tube  radius  in  Region  B  Is  demon- 
strated in  Fig.  6.5.   The  drift-tube  radius  in  Region  B  should  be  increased 
because  the  best  amplification  is  achieved  for  smaller  values  of  b/a.   It 
is  apparent  from  Fig.  6.5  that  the  decrease  in  reduction  factor  resulting 
from  the  increase  in  drift-tube  potential  can  be  compensated  at  least 
partially  by  increasing  the  drift-tube  radius.   Amplification  is  achieved 
in  the  third  section  if  the  drift-tube  potential  is  decreased  at  the  sec- 
ond transition.   The  amplification  can  be  further  enhanced  also  by  restor- 
ing the  drift-tube  radius  to  the  value  of  Region  A. 

The  principle  purpose  of  the  first  transition  is  to  change  the  drift- 
tube  potential  from  a  value  Va  to  a  higher  value  V^,  without  decreasing 
the  velocity  amplification,  or  perhaps  to  even  increase  It  somewhat  by 
increasing  the  drift-tube  radius.   The  drift-tube  potential  is  returned 
to  the  original  value  (Va)  at  the  second  transition.   The  amplification 
resulting  from  the  decrease  In  drift-tube  potential  was  the  subject  of  a 
preceding  section. 

The  locations  of  the  two  transitions  are  quite  important.  Small- 
signal  theory  predicts  that  each  drift  tube  should  be  a  quarter  plasma 
wavelength  long.   The  velocity  at  the  second  transition  is  shown  in 


3h 


05  - 


,ok    ■■ 


03  - 


,02   - 


.01   ■ 


Fig.  6.5  -  Velocity  In  Region  B  as  a  function  of  the  drift-tube  radius, 
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Fig.  6.6  for  the  case  where  the  velocity  is  maximum,  i.e.,  the  transition 
locations  are  optimized,  and  for  the  case  where  each  drift-tube  section 
is  a  quarter  plasma  wavelength  long.   The  maximum  velocity  obtained  In 
Region  A  is  also  plotted  for  comparison.   It  is  obvious  that  the  velo- 
city modulation  can  be  improved  by  optimizing  the  transition  locations. 
The  optimum  locations  of  the  first  and  second  transitions  as  a  function 
of  cc   are  shown  in  Fig.  6.7  and  6.8  respectively. 

If  the  second  transition  is  located  at  the  point  of  maximum  a-c 
velocity  in  Region  B,  then  the  maximum  velocity  in  the  third  region  is 


obtained  by  multiplying  the  velocity  at  the  transition  by  \l  V.  /V  .   The 
ratio  of  the  maximum  velocity  in  Region  C  to  the  maximum  velocity  in 
Region  A  is  plotted  in  Fig   6.9  for  the  solutions  predicted  by  the  small- 
signal  theory  and  by  the  nonlinear  theory  for  the  two  aforementioned 
cases.   The  nonlinear  theory  definitely  shows  the  effect  of  saturation. 
Although  it  is  possible  to  improve  the  gain  by  optimizing  the  drift-tube 
lengths,  it  is  apparent  that  for  practical  values  of  a   the  amplifier  is 
useless.   Perhaps  for  larger  voltage  jumps,  the  amplifier  will  not 
saturate  at  such  a  low  value  of  a.      Such  a  conclusion  can  only  be  veri- 
fied by  further  computer  work,  and  of  course  in  the  final   analysis,  by 
experimental  Investigation. 
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Fig.  6.6  -  Maximum  velocity  in  Region  B  as  a  function  of  a. 
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(b/a)B  =  1/2 


0.2 


0.3    a 

Fig.  6.7  -  Location  of  the  first  transition  for  maximum 
velocity  in  Region  B. 


VVa  =  2-0 
(b/a)B  =  1/2 


0.2 


0.3    a 


Fig.  6.8  -  Location  of  the  maximum  velocity  in  Region  B. 
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Smal 1 -s  ignal 
predict  ion 


Each  transition 
made  at  optimum 
pos  i  t  ion 


Each  transition 
at  a  quarter 
reduced  plasna 
wavel ength 


Fig.  6.9  -  Effect  of  drift-tube  lengths  upon  the  velocity 
ampl i  f  i  cat  ion. 


CHAPTER  7 
SUMMARY  AND  CONCLUSIONS 

Tien  and  Field;  Birdsall  and  Whinnery;  Peter,  Bloom  and  Ruetz;  and 
others  had  utilized  smal 1 -s ignal  analyses  to  investigate  the  performance 
of  a  space-charge-wave  amplifier  which  employed  changes  in  the  potential 
and  radius  of  the  drift  tube  at  the  paint  of  minimum  and  maximum  velo- 
city variations  to  enhance  the  amplification.   The  amplifier  was  found 
to  have  excellent  amplification  properties  which  depended  only  upon  the 
ratio  of  the  potential  change  between  sections,  and  to  a  lesser  extent, 
upon  the  change  in  drift-tube  radius.   Experimental  investigation,  how- 
ever, had  shown  that  the  amplifier  definitely  saturated  at  quite  low 
signal  levels,  and  therefore  exhibited  low  efficiency. 

The  purpose  of  this  dissertation  was  to  study  the  saturation  pheno- 
mena, and  try  to  determine  if  the  performance  could  be  improved  by  ad- 
justing the  lengths  of  the  drift-tube  sections.   A  nonlinear  theory 
proposed  by  Paschke  for  the  study  of  the  nonlinear  effects  in  a  basic 
two-cavity  klystron  was  utilized.   The  nonlinear  space-charge-wave  eq- 
uation was  developed  in  Chapter  3,  and  the  general  solutions  were  ob- 
tained in  Chapter  k   by  employing  the  method  of  successive  approximations 
as  proposed  by  Paschke.   In  the  following  two  chapters  particular  solu- 
tions for  the  first  and  second  drift-tube  sections  were  obtained  by 
applying  boundary  conditions  of  conservation  of  energy  and  continuity 
of  current  to  the  transition  locations. 
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The  solutions  for  the  first  region  applied  to  the  conventional  two- 
cavity  klystron.   Solutions  for  current  and  velocity  were  plotted  in 
Fig.  5-1  and  Fig.  5-2  for  three  values  of  the  voltage  modulation  index 
(a).   These  solutions  obtained  from  the  nonlinear  theory  demonstrated 
the  saturation  phenomena,  and  also  showed  that  the  locations  of  minimum 
and  maximum  current  and  velocity  deviated  significantly  from  the  multiples 
of  quarter  plasma  wavelengths  predicted  by  the  small-signal  theory.   It 
was  shown  in  Fig.  5-^  that  the  locations  of  the  minimum  and  maximum  cur- 
rent for  a   =  0.3  are  only  approximately  two-thirds  of  the  distances  pre- 
dicted by  small-signal  theory.   Consequently,  if  the  first  transition 
were  placed  a  quarter  plasma  wavelength  from  the  buncher  it  would  not  be 
at  the  position  of  minimum  velocity,  nor  would  a  transition  placed  one- 
half  of  a  plasma  wavelength  be  at  the  point  of  maximum  velocity  variation. 

The  solutions  for  the  second  region  gave  information  concerning  the 
current  and  velocity  after  the  space- charge-wave  had  undergone  a  velocity- 
jump.   Several  different  combinations  of  voltage  ratios,  drift-tube  radii, 
and  transition  locations  were  explored.   These  studies  fell  into  two 
general  catagories;  a)  transition  near  a  velocity  maximum,  and  b)  transi- 
tion near  a  velocity  minimum. 

For  a  transition  near  a  velocity  maximum  it  was  determined  that  the 
best  amplification  was  obtained  if  the  voltage  and  radius  of  the  second 
drift  tube  were  decreased.   The  major  effect  was  due  to  the  v  1 tage  tran- 
sition rather  than  the  change  in  the  drift-tube  radius,  as  can  be  seen 
by  comparison  of  Fig.  6.1  and  Fig.  6.2.   It  was  demonstrated  in  Fig.  6.3 
that  the  best  position  for  locating  the  transition  was  not  at  the  point 
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of  maximum  velocity  variation,  but  s  me  twenty-two  degrees  closer  to  the 
input.   For  larger  values  of  a.   the  shift  would  most  likely  be  even  greater 
since  the  bean  acts  in  a  more  nonlinear  fashion. 

It  was  also  demonstrated  in  Fig.  6.1  that  if  the  amplifier  utilized 
only  two  drift-tube  sections,  the  current  output  was  not  as  good  as  a 
conventional  two-cavity  klystron  which  had  the  same  signal  voltage  applied 
to  the  input,  but  utilized  the  lower  drift-tube  voltage,  and  consequently 
the  higher  beam  perveance,  of  the  second  section.   Consequently,  it  was 
necessary  to  utilize  an  increase  in  potential  at  the  first  transition  in 
order  to  insure  that  the  same  maximum  beam  perveance  was  available  for 
the  first  section  of  both  amplifiers.   It  was  demonstrated  in  Fig.  6.4 
and  Fig.  6.5  that  a  single  transition  near  a  velocity  minimum  with  an  in- 
crease in  potential  did  not  produce  a  useful  amplifier  either.   However, 
combining  this  transition  with  a  second  transition  near  a  velocity  maximum, 
wi  ere  the  drift-tube  potential  was  decreased  to  the  original  value,  pro- 
duced an  amplifier  from  which  gain  Could  be  obtained. 

A  comparison  was  made  in  Fig.  6.6  of  the  velocity  amplitude  for  a 
conventional  klystrn,  and  for  the  space-charge-wave  amplifier  first  using 
quarter  plasma  wavelength  sections,  and  then  for  the  case  where  the  drift- 
tube  lengths  had  been  optimized.   It  was  evident  that  the  amplifier  had 
improved  performance  for  larger  input  signals  only  if  the  drift-tube 
lengths  were  optimized.   The  required  transition  locations  were  shown  in 
Fig.  6.7  for  the  first  transition,  and  in  Fig.  6.8  for  the  second  transition, 
It  was  apparent  that  these  locations  were  definite  functions  of  a,    and  that 
the  locations  shifted  toward  the  input  as  ct    increased.   Comparison  of  Fig. 
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6.7  and  Fig.  5-^  indicates  that  tne  first  transition  location  definitely 
does  not  coincide  with  the  point  of  minimum  velocity  variation  as  predicted 
by  the  small -signal  theory,  and  in  fact,  closely  follows  the  location  of 
the  maximum  current. 

The  ratio  of  maximum  velocity  in  the  third  region  to  the  velocity 
initially  introduced  by  the  input  signal  was  plotted  in  Fig.  6.9  for  the 
small -signal  theory  and  the  nonlinear  prediction  for  drift  tubes  with 
quarter  plasma  wavelength  sections  and  optimized  sections.   It  is  apparent 
that  the  nonlinear  theory  indicates  the  saturation  effect  which  was  not 
predicted  by  small-signal  theory,  and  that  performance  can  be  improved 
by  optimizing  the  lengths  of  t:ie  drift  tubes.   However,  it  is  also  apparent 
that  the  improvement  is  not  adequate,  since  the  ratio  still  falls  below 
unity  for  reasonably  snail  input  signals. 

It  should  be  pointed  out  again,  however,  that  the  theory  utilized  in 
this  analysis  is  valid  only  to  the  point  of  electron  crossover.   For  the 
parameters  utilized  in  this  invest igat ion  such  crossover  takes  place  for 
a   =  0.15,  i.e.,  v/el  1  below  the  point  where  saturation  effects  are   a  serious 
problem.   Consequently,  the  saturation  phenomena  could  not  be  c ■mpletely 
analyzed  by  this  theory.   It  is  also  possible  that  gain  could  have  been 
achieved  for  larger  values  of  a  by  utilizing  a  larger  voltage  jump.   The 
two  to  one  change  in  potential  used  for  several  of  the  curves  is  well  below 
the  changes  utilized  by  Tien  and  Field  in  their  original  work.   Only  further 
computer  investigations  for  various  possible  parameter  combinations  and 
detailed  experimental  work  can  yield  a  more  complete  picture  of  the  satura- 
tion phenomena. 


APPENDIX  A 

LIST  OF  PRINCIPLE  SYMBOLS 

a    Drift-tube  radius,  and  subscript  used  to  indicate  a  variable  in 
Region  A 

b    Beam  radius,  and  subscript  used  to  indicate  a  variable  in  Region  B 

-19 
e    Electron  charge  (-  1.6  x  10   coulomb) 

f    Electric  flux  coupling  beam  and  drift  tube 

Fundamental  component  of  current  density 

Electron  current  density  in  unmodulated  beam 

Electron  current  density  in  modulated  beam 

A-c  electron  current  density  (it  -  iQ) 

Normalized  electron  current  density  (i/i0) 

m    Electron  rest  mass  (9-1  x  10"^  ki logram) 

n     Integer  used  to  represent  the  harmonic  number 

s    V6q2  =  Rl/R2 
t    Time 


Electron  velocity  in  unmodulated  beam  -J 2rjV0 

Electron  velocity  In  modulated  beam 

A-c  component  of  electron  velocity  in  modulated  beam  (v^  -  vQ) 

Normalized  component  of  a-c  velocity  (v/vQ) 

Axial  distance  measured  from  input  cavity 

Normalized  axial  distance  (B  ,  z) 

v  qla  ' 
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z    Component  of  electron  polarization  obtained  from  n   order  solution 

z    Polarization,  or  displacement,  of  an  electron  from  its  position  in 
the  unmodulated  beam 

A    Beam  cross  sectional  area 

E    Electric  field  resulting  from  space  charge 

J,    Bessel  function  of  1st  order  and  1st  kind 

6   3/2 
K    Beam  micro-perveance  ( lQ  10  /V    ) 

M    Displacement  current/unit  length  coupling  drift  tube  and  electron 
beam 

R  Plasma  reduction  factor  at  frequency  noo 

V  Amplitude  of  a-c  potential  applied  to  input  cavity 

V  Drift-tube  potential  in  Region  N 
X  Ballistic  bunching  parameter 

X    Bunching  parameter  including  space-charge  effects 

a.  Voltage  modulation  index  (Vm/Voa) 

S    Electronic  wave  number  (uj/v  ) 
e  v    o' 

B  b  Normalized  beam  radius 
e 

Sp   Plasma  wave  number  (<*>D/v0) 

8    Reduced  plasma  wave  number  (Rn8n) 

qn  "  k 

-12 
€    Permittivity  of  free  space  (8.85^  x  10    farads/meter) 

p  Electron  charge  density  in  unmodulated  beam 

pt  Electron  charge  density  in  modulated  beam 

p  A-c  component  of  charge  density  in  modulated  beam 

r\  Ratio  of  lel/m  (1.76  x  10   coulomb/kilogram) 

w  Angular  frequency  of  signal  voltage  to  Input  cavity 

u  Angular  frequency  of  plasma  oscillation  for  a  beam  of  infinite  size 

(mpJ/0,/2 
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oj    Reduced  plasma  frequency  of  oscillation  (R|to  ) 

\x  Gap  coupling  coefficient 

0^        Phase  shift  of  fundamental  component  of  current  density 


APPENDIX  B 

LIST  OF  PARAMETERS  USED  IN  EQUATIONS 

The  following  group  of  parameters  is  needed  in  order  to  express 
the  parameters  introduced  in  (4.11)  in  a  compact  fashion  for  computer 
analys  i  s . 

x   =  3e/^ql 


s 

"    V  V 

t 

B    V9q2 

g 

=   (1    -  4s2)"1 

bi 

=    1    -   s 

b2 

=   1   +  s 

b3 

=  2s   -    1 

bk 

=   2s  +   1 

BOJ 

"    3ql/3qla 

BQ2 

-   aq2/3qla 

BE 

"    0e/0qla 

Pl 

=   A   B  +   C   D 

P2 

=   A2   +   B2   +   C2   +   D2 

P3 

=  A   C  +   B   B 

p4 

-2        -2        -2 

-  A    -  b    +  £    - 

-D2 

p5 

-  a  6  +  §  e 

Pa 

=  A2  +   §2   -   C2    • 

-   B2 

106 


107 


P7  =  A  B  -  C  5 

P8  =  A2  -  B2  -  C2  +  D2 

P9  =  A   C  -  6  6 

P,«  =  A  5  -  B  C 


The  following  parameters  are  introduced  in  (4.11).   The  bar  (  ) 
indicates  the  parameter  Is  normalized  by  dividing  by  (3Q|  . 
K,   =1-5  (BOJ)2  P, 
R2  =  -  0.75  (BQI)2  P^ 
R3  =  x  (BQI)2  P3/2 
Kj,  =1-5  (BQI)2  (P7  -  x  Pg) 
K5  =  -  0.75  (BQI)2  (P8  +  2  x  P5) 
R$  =  -  0.25  (BQI)2  x  P6 
K?  =  0.75  (BQI)2  (x  P8  +  2  P5) 
K8  =1.5  (BQI)2  (x  P?  -  Pg) 

The  following  parameters  are  needed  for  describing  the  parameters 
introduced  in  (4.20). 
X,   =  2x2  +  3 
X2  =  x2  +  1 
X7  =  2x2  -  1 


X,   =  2x2  -  7 


,2 


X6  =  x2  -  1.5 
X?  =  x2  -  3 
T,   =  kt2   +  1 
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T2 

=  4tz  -  1 

Ha 

=  BQI  (X2  A  +  2  x  D) 

Hb 

=  BQI  (X2  B  -  2  x  C) 

Hc 

=  BOJ  (X2  C  -  2  x  B) 

Hd 

=  BQI  (X2  6  +  2  x  A) 

He 

=  kt2   +  x2  -  6s2  -  3 

H* 

=  kt2  +   x2/2  -  3 

H:        = 


=  6s2g  X,    +  x2   X2 


=   s2   g  X,    +  x2/3 

s2  g  X/,  +  x2/3 
=  g   (8t2   -  3)/2  +   1/2 


c2/6 


c 
Kd 
Ke 
Kf 

Kj 

KJ 
Kk 

•Sn 

c2 


=  BQI 

=  BQI 

=  BQI 

=  BQI 

=  BQ2 

=  BQ2 

=  BQ2 

=  BQ2 

=  BQ2 

=  BQ2 

-  BQ2 
=  BQ2 

-  BQI 

-  BQI 


x  A  +   D) 
x   B   -   C) 

x  C  -  B) 

x   D   +  A) 

2    t    T    -  M) 

2   t   M   -    T) 

2   t    L   +   J) 

2   t   J  +   [) 

k  x   T   -  Tj    M/s) 

k  x   J  +  L  Tj/s) 

-k  x  M  +  T  Tj/s) 

4  x   L  +  J  Tj/s) 

X6  R2  +  Hg   K5/6  +  s    t   g  X4   K7  +  x   (BE)2  X5   P5A)/2 

X6  R,    +  Hh   Kk  -   s    t   g  X4  R8  +   (BE)2  X5  P?A)/2 
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C3 
c5 

C7 

c8 

Cn 


10 


1  1 


12 


13 

14 
;15 
;16 


17 


19 


'20 


BOJ    (2  x  X7   K|/3   +  x  Hj    K4  +  s2   g  X]    Kg  +   (BE)2  X5  PgA)/2 
-BQ.1    (2   x  X7   R2/3   +  x   H?    K5    -   s2   g   X]    R7  +   (BE)2   X5   P5A)/2 

-  (x   R2  +  2    s    t    g   R5  +  Hj    R7)/2 

-  (x   Rj    +  2   s   t   g   R/+   -  Hj    K8)/2 

-  (X?   K,/3    -   Hj    Kk  +  2   s    t   g   Rg)/2 
(X7  K2/3   -  H.   Kg  -  2  s  t  g  K?)/2 

g   s    t    (R7   Ha    -    Rg   Hb   -   R^   Hc   +   R5   Hd) 

-  s    t   g    (R5   Ha   +   R^   Hb   +   Kg   Hc   +   R?   Hd) 
BOJ    (x   Kj   +  1.5  X3    Kf)/2 

BQ.1    (-  x   K.    +   1  .5  X3    K  )/2 

BOJ    (x   K      -    1.5   X,    K.  )/2 
m  3      n 

BQ1    (x   Kk   -    1.5  X3    Ke)/2 

(x   Kf   -    1.5   BQ2  T2  M/s)/2 

(x   K     -    1.5   B0.2  T2    L/s)/2 

(-  x   Kh  +   1.5   B0.2  T2   J/s)/2 

(-  x   Ke  +   1.5   BQ2  T2   T/s)/2 

s    t   g    (R7   Ha   +   Rg   Hb   -    Kk   Hc   -    R,-   H<J) 

s    t   g    (Rg   Ha    -   R?  Hb   -    K5   Hc   +   Kk  H,) 

s    t   g    (R4  Ha  +   R5   Hb  +  R7  Hc  ♦   Rg  Hd) 


C22  =   s    t   g   (R5   Ha   -   Kk  Hb  +  Rg  Hc   -   K?  Hd) 
=   5    t   g    (K8   Ha   +   R7  Hb   -   K5   Hc   -   K4  Hd) 
=   C3    -   2   BE  X      Rj/3 
C^  +  2    BE   X7   K^/3 

-  C,    +   BQ1    X6   R2 

-  C2  +  BOJ    X$   K] 


C23 

C25 
C26 

C27 
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C28  =  C7  +  X7  R,/3 

^29  =  ^8  "  ^7  ^2 

C30  =  "  C5  "  x  K2 
C3,   =   -  C6  -  x  R, 

C32  =   s    t   g    (-   R4  Ha  +  K5   Hb   -   K7  Hc  +  K8  Hd) 

The   following   parameters   are    introduced    in    (4.20).      The   bar    (    ) 
indicates   that    the  parameter    is   normalized   by   dividing  by   (3  j    . 
Kg     ■   BE(Kg  A   -   K3    C)He/2   +   («c   i<3    -    Ka    Kg)    Hf   -   2    BQJ    Xg  A   E   +   2   x   Ka   E 

+  A  C|    -   B   C2   +  C   C3    -   6   C^  +   Ka   C5   -   Kb   Cg  +  Kc   C7   -   Kd   C8  +  C]9 

+  2   s    t    (Kg   Ha    -    K3    Hc) 
K,0  =  A  Cn    +  B   C]2  +  C   C,3      +  6   C]Z+  +  Ka   C,5  +  Kb   C]6  +  Kc  C]?+  Krf  C,g 

+  BE   (Ha  M  +  Hb    L   -   Hc   J   -   Hd   I) 
K,,    =  A  C,,    -   B   C12  +  C   C,3    '   0   C]k+   Ka   C,5   -   Kb  C|6  +  K£   C]?   -   Kd   C)Q 

+   BE    (Ha   M   -   Hb    L   -   Hc   J   +   Hd    I) 
K)2  =  A   C]    +  B   C2  +  C   C3   +  6  C4  +  Ka   C5  +   Kb   Cg  +  Kc   C7  +  Kd   C3  +  Cg 
K]3   =    BQ1       (x   He   Kg/2      -   2   Xg   E)    B   +   K?    (Kd   Hf    -   BQJ    x   Hg   5/2) 

+   Kb   (2   x   E   -   Hf   Kg)    +   C20   +  2   s    t    (Kg   Hb   -    K3    Hd) 

-  (A  C2  +  B  Ci  +  C  C^  +  5  C3  +  Ka  Cg"  +  Kb  C5  +  Kc  C8  +  Kd  C7) 
K|i,  =  AC,2-BC,i+C  C14  -  6  C)3  +  Ka  C]6  "  Kb  C]5  +  Kc  C]8  -  Kd  C]? 

+  BE  (Ha  L  -  Hb  M  -  Hc  T  +  Hd  J) 

K)5  =  A  C]2  +  B  C,,  +  C  C]k   +  6  C,3  +  Ka  C,g  +  Kb  C,5  +  Kc  C,g  +  Kd  C,? 

+  BE  (Ha  L  +  Hb  M  -  Hc  T  -  Hd  J) 

K|6  ="A  C2  +  B  Ci  -  C  C4  +  5  C3  -  Ka  C6  +  Kb  C5  -  Kc  C8  +  Kd  C;  +  C23 

K,7  =  A  C2k   -    B  C25  +  C  C2g  -  6  C2?  +  Ka  C28  -  Kb  C2g  +  Kc  C3Q 

-  Kd  C31  -  2  BQJ  Xg  C  E  -  (1/2)BE  (A  K3  +  C  Kg)  Hg  +  (Ka  K^ 


Ill 


+  Kc  i^)  Hf  +  2  x  Kc  E  -  2  s  t  (K3  Ha  +  R6  Hc)  -  C2, 
Kl8  =  A  C,3  +  B  C]k   -  C  C,,  -  5  C]2  +  Ka  C]?  +  Kb  C]g  -  Kc  C,5  -  Kd  C]6 

-  BE  (Ha  J  +  Hb  T  +  Hc  M  +  Hd  L) 

K,g  =   K]8  +  2    (-   B   C1Z+  +   6  C,2    -   Kb   C,8  +  Kd   C,6)    +  2    (Hfa    I    +  Hd   [)    BE 
K20  =   A   C2k  +   B   C25   +   C   C26  +   5   C27  +   Ka   C2g   +   Kb   C2g   +   Kc   C30   +   Krf   Cj, 

+  c32 

Rj,  -  -  (A  C25  +  B  C2k  *   C  C2?  +  5  C26  +  Ka  C29  +  Kb  C28  +  Kc  Cj, 

+  Kd  C30)  -  BQ1  2  X6  5  E  -(1/2)  BE(K3  B  +  5  K6)  He  +  (R3  Kb 

+  K6  Kd)  Hf  +  2  x  Kd  E  -  2  s  t  (K3  Hb  +  i^  Hd)  -  C22 

K22  =  A  C]k  -   B  C,3  -  C  C,2  +DC,,  +  Ka  Cl8  -  Kb  C1?  -  Kc  C,6  +  Kd  C,5 

-  BE  (Ha  T  -  Hb  J  +  Hc  L  -  Hd  ft) 

K23   =   K22   +  2    (g   C|3   -   6  CH   +   Kb   C,7   -    Kd   C,5)    -   2    BE    (Hb   J   +   Hd   M) 
K2l^  =    -  A   C25   +   B   C2Z+   -    C   C27  +   6   C26   -    Ka   C2g   +  Kb   C28   -   «c   C3 , 

+   Kd   C30   +   C10 

The  following  functions  of  z_  are  needed  in  (1+.28). 
Fa(z)  =  ^(z)  I  B  cos  (3q,z)  -  A  sin  (  Jqlz)]  +  f2(z)  [5  cos  (3q]z) 

-  C  sin  (3q]z)] 

Fb(z)  =  f2(z)  [B  cos  (3q]z)  -  A  sin  (3q]z)]  +  f?(z)  [5  cos  (3q,z) 

-  C  sin  (3a]z)] 

Fc(z)  =  -  f4(z)  [Kc  cos  (3q,z)  +  Kd  sin  (3qlz)]/2  +  f  (z)  [Ka  cos(3qlz) 

+  Kb  sin  (3q)z)]/2 
Fd(z)  =  -  f5(z)  [Kc  cos  (3q,z)  +  Kd  sin  (3q,z)]/2  -  f^z)  [Ka  cos  (Sq]z) 

+  Kb  sin  (3q]z)]/2 
where 
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f](z)  =  BQ1  [2E  -  2  s  t  Kg  +  {2  s  g  (K5  -  x  K?)  -K2)  cos  (2pq]z) 
+  {K,  -  2  s2  g  (R4  +  x  Kg)}  sin  (2Sq)z)  -  Kf  cos  (3q2z) 

-  Kg  sin  (3q2z)]/2 

f2(z)    =   BQ1    [2   s    t    K3   +  2   s2   g    ((x   K^  +  Kg)    cos    (23q]z)    +   (x   i<5   -   K?) 

.sin    (23q,z)}    +  Kh   cos    (3q2z)    +  Ke  sin    (3q2z)]/2 
f3(z)   =    BQ1    [2  E   -  R2  COS    (23q)z)    +  K,    sin    (23q,z)]    -   f,(z) 
f4(z)    =   2   E  +  BE   BQ1    P,0  +  (i'2)(BQl)2   ?2   -   l^/x 

+  (1/3)    [g    (1    +  2   s2)    R5   -   K2]    cos    (26qlz) 

-  (1/3)  [g  (1  +  2  s2)  Kk   -   R,]  sin  (23q]Z) 
+  BQ2  [J  cos  (3 _2z)  -  J  sin  (3 «z) ] 

f5(z)  =  K3/x  +  (1/3)  g  (1  +  2  s2)  [Kg  cos  (23q]z)  -  Ky  sin  (23q]z)] 

+  BQ2  [L  cos  (3q2z)  -  M  sin  (3q2z)] 
f6(z)  =  ^(z)  -  4  E  -  2  BE  BQ.1  P]Q  -  (BQJ)2  P£ 

+  2  [R2  cos  (23  ,z)  -  K,  sin  (23  ,z)]/3 

The   following   constants   are  associated  with   the   third-order   solu- 
tions   in   Region  A. 
A,      =   -  A2   -  A3    -  A^ 

A2     =(3/128)    ga   ta[8   s/    (b2a   -   s/)    -(6sa+,)]/b3a 
A3     =   (3/128)   ga   ta    [8  s/   (b,a   -  sa2)   +  (6  sa   -   1)]   /b^ 
A4     =   -    (15/256)    sa   ta  ga 
A5     =   -   (12/5)    A4 
A6     »    15/64  -  A7  b,a/sa   -  Ag   b2a/sa   -  3   Ag   -  A,0 

A7     =   "   Sa  bla  %    [9   <•    "  2   Sa^   +  k  Xa   'a   b2a   <b4a   "  2   Sa^128   b3a> 
A8     =   sa   b2a  9a    ^    ('    "  2   sa2)    +  4  *a   <a   bla    <b3a  +  2   -a^^^28   bO 
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Ag 

A10 

All 


-  ga    [9    (10   sa2    -l)+2  xa2    (7   sa2    -    l)]/102if 

-  9a    tt   2    (7   -  8   s   2)    +  27   s   2  +  x  2   -  9/2J/128 

a  o  a  a  a 


=  Ac 


A12  = 


'13 


[)k 


M5 


A5   -  A,    +  xa   ga    (1    -  28   s/)/128 

*a  9a    t2  sa   (7  -    10   sg2)    +  3    (I    -  2   sa2)]/128     -  A2   b|a/sa 
xa  ga   [2sa  (7  -   10  sa2)   -3(1-2  sa2) J/128     -  A3  b2a/sa 
xa  ga  (1   -  16  sa2)/128  -  3  A4 


A16="A 
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>17  =  A6  +  A]0  +  ga   (1    +  8  sa2)   X2/128  -   (5   -  2   ta2)/64 

'18  =  A7b.a/Sa   "  *a    [l  *a   la   < '    "   Sa^    \a  +  3    <'    "  2   Sa^  ]/64 
Aftbo/S;,    -a      [2x     t      (I-   s  a2)    b,      +  3    ( 1    -  2   s   2)  ]/6k 


k19 


*20 


A8u2a/ba 

3  Ag   -    (3   +  xa2   -    12  ga   sa2)/128 


3a 


APPENDIX  C 
DISCUSSION  OF  ELECTRON  OVERTAKING 

The  space-charge-wave  equations  derived  from  the  plasma  model  are 
valid  only  so  long  as  the  electron  velocity  is  a  single-valued  function 
of  position,  consistent  with  the  polarization  system  of  coordinates. 
The  problem  remains  to  identify  just  when  this  limitation  is  violated. 
Physically,  the  location  at  which  overtaking  occurs  is  the  point  where 
the  change  in  polarization  (z  )  exceeds  the  change  in  position  along  the 
beam  axis  (z).   Hence,  overtaking  occurs  when  the  slope  of  zp  verses  z 
is  greater  than  one. 

Engler  [26]  has  investigated  the  problem  on  the  assumption  that 
the  polarization  can  be  represented  by  the  first-order  solution  (z,) 
alone.   This  solution  is  obtained  by  substituting  (5-15)  and  (5-17)  into 

(*».5). 

5]  ■  (a/2)    sin(z)  sin(wt  -  xa  z)  (C.l) 

Taking  the  partial  derivative  of  (C.l)  with  respect  to  I  gives 

i|   =  (a/2)    [  -  xa  sin(z)  cos(wt  -  xa  z) 

+  cos(z)  sin(ut  -  xa  5)3.  (C.2) 

Equation  (C.2)  may  be  written  as  a  single  function  of  time  in  the  form 
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i|  '  =  (a/2)>J  x2   sin2(I)  +  cos2(i)   cos(wt  -  x  I  -  <t>) .      (C3) 

_  i 
Overtaking  first  occurs  when  Zi   =1.   Since  xa  >  1 ,  the  coefficient  of 

cos(u)t)  can  be  maximized  for  a  given  value  of  a  if  z  =  90°.   Thus  the 

maximum  value  of  a    Is 


"max  "  2/V  <C^ 


Engler  mentions  that  if  higher-degree  terms  for  the  polarization  are 
retained,  the  value  of  a  may  be  raised. 

To  confirm  this  conclusion  all  three  orders  of  terms  in  the  polari- 
zation expression  are  retained.   The  total  polarization  is  the  sum  of 
C+-5),  (^.17)i  and  (4.21),  where  the  results  of  Chapter  5  are  utilized 
to  evaluate  the  coefficients.   Since  the  actual  polarization  can  be 
found  only  for  explicit  values  of  time,  the  polarization  is  computed  for 
cot  =  0°  and  u>t  =  90°. 

Fig.  C.l  shows  how  the  polarization  varies  as  a  function  of  a. 
The  saturation  effect,  which  is  present  for  the  current  envelope,  is 
also  demonstrated  here.   The  important  point  to  note  is  the  fact  that 
the  magnitude  of  the  polarization  and  the  portion  of  the  cycle  for  which 
overtaking  occurs  both  rise  with  increasing  values  of  a.      Fig.  C.l 
indicates  that  overtaking  is  present  over  most  of  the  cycle  for  a   =  0.3. 

Fig.  C.2  shows  the  change  in  polarization  with  time.   Although  the 
polarization  at  any  point  on  the  beam  definitely  changes  with  time,  the 
overall  pattern  seems  to  be  nearly  the  same  regardless  of  the  value  of 
time  chosen.   The  patterns  are  not  identical,  however,  as  the  second- 


116 


Fig.  CI  -  Polarization  as  a  function  of  <£. 


a  =  0.15 
K  -  1.0 

b/a  »  2/3 
Seab=  1-° 


Fig.  C.2  -  Polarization  as  a  function  of  ut, 
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order  term  (z~)  contains  a  t ime- independent  function  of  position.   Thus 

It  may  be  concluded  that  overtaking  will  be  spotted  for  a  given  value 

of  a   regardless  of  the  time  value  used;  however,  the  location  along 

the  beam  at  which  overtaking  occurs  is  a  function  of  the  time. 

Careful  examination  of  Fig.  C.2  shows  that  overtaking  occurs  In  the 

vicinity  of  92°  for  cot  =  0.   With  the  set  of  parameters  used  in  Region  A, 

(C.4)  predicts  that  amax  =  0.175-   Hence,  by  including  higher-order  terms, 

overtaking  is  shown  to  occur  for  a  smaller  value  of  a  rather  than  a 

max 

larger  value  as  predicted  by  Engler.   It  is  also  apparent  that  the  degree 
of  overtaking  is  slight,  and  the  value  of  o^gx  predicted  by  (C.4)  is  an 
acceptable  guide  to  follow. 


APPENDIX  D 

METHOD  FOR  OBTAINING  SOLUTIONS 
FOR  THE  CONSTANTS  OF  INTEGRATION  IN  REGION  B 

In  this  appendix  the  procedure  for  solving  four  simultaneous 
equations  using  Cramer's  rule  is  discussed.   In  Chapter  6  solutions 
for  the  first-,  second-,  and  third-order  constants  of  integration, 
which  are  required  to  obtain  particular  solutions  in  Region  B,  are 
each  left  in  the  form  of  four  simultaneous  equations.   The  first-order 
constants  are  found  in  (6.14)  and  (6.I7),  the  second-order  in  (6.24) 
and  (6.31),  and  the  third-order  in  (6.39)  and  (6.44). 

Consider  the  following  system  of  equations: 

Fj  =  C,  S,  +  C2  S2  +  C3  S3  +  C^  S^  (D.la) 

F2  =  C5  Sl  +  C6  S2  +  C7  S3  +  C8  S4  (D.lb) 

F-,  "  c«  s,  +  c,«  So  +  c,,  s-,  +  C,0  S.  (D.lc) 

3    9  1    10  2    11  3    12  4  v 

F4  =  C13  Sl  +  C14  S2  +  C15  S3  +  C16  V        (D,,d) 
The  solutions  to  this  set  of  equations  are 

Sl  =  [F1  Mll  "  F2  M21  +  F3  M31  ■  F4  M41]/D        (D'2a) 
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52  =  [  -  F,  M|2  +  F2  M22  -  F3  M32  +  F^  M^2]/D      (D.2b) 

53  =  [F,  M,3  -  F2  M23  +  F3  M33  -  F^  M^l/D         (D.2c) 

Sk  =    [    ~    Fl  MI4  +  F2  M24  "  F3  V  +  Fk   MW/D-      (D'2d) 

The  symbol  D  represents  the  fourth-order  determinant  formed  by  the 
coefficients  C]  through  C]£.   The  symbol  M^  represents  the  minor  formed 
when  the  mth  row  and  n_th  column  are  deleted  from  D.   In  general   these 
solutions  would  appear  to  be  quite  complicated,  since  there  are  sixteen 
unique  minors  which  must  be  evaluated.   In  addition  there  are  three 
sets  of  such  equations  which  must  be  evaluated,  one  for  each  order 
solution  in  Region  B. 

The  work  is  considerably  simplified,  however,  when  the  equations 
which  (D.l)  represent  are  examined.   The  first  important  point  is  the 
fact  that  for  a  particular  order  solution  there  are  only  four  unique 
constants.   These  constants  are  the  products  of  sine  and  cosine  terms. 
For  the  first-order  solution  they  are 

C,  =  cos(BQlb  Za)  sin(BEb  Za)  (D.3a) 

C2  =  sin(BQlb  Za)  sin(BEb  Za)  (D.3b) 

C  =  cos(BQJb  Zg)  cos(BEb  Za)  (D.3c) 

C^  =  sin(BQlb  Za)  cos(BEb  Za) .  (D.3d) 
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The  four  first-order  equations  may  be  put  in  the  form 


F,  =   C]  Sj  +  C2  S2  +  C3  S3  +  C^  S^  (D.^a) 

F2  =  -  C2  S,  +  C,  S2  -  C/^  S3  +  C3  S^  (D.4b) 


F3  =  -  C3  S,  -  C^  S2  +  C,  S3  +  C2  S^  (D.4c) 


Fk   =   Ck   Si  "  C3  S2  "  C2  S3  +  Cl  V  (D'^d) 


The  solutions  S.,  S  ,  S  ,  and  S,  refer  to  the  first-order  constants 
A, ,  B. ,  C. ,  and  5,  respectively.   The  forcing  functions  are   given  by 

F,  =  (a/2)  (xa/BEb)  sin(Za)  cos(xa  Zg)  (D.5a) 

F2  =  (a/2)  (Vg/Vb)  ( l/BQlb)  cos(Za)  cos(xa  Za)  (D.5b) 

F3  =  (a/2)  (xa/BEb)  sin(Za)  sin(xa  Za)  (D.5c) 

F^  =  (a/2)  (Va/Vb)  (1  BQlb)  cos(Za)  sin(xa  Za) .         (D.5d) 

The  solutions  to  (D.k)    are 

S,  =  [F,  M,,  +  F2  M,2  -  F3  M,3  -  F^  M]i+]/D  (D.6a) 

S,  =  [-  F,  M,,  +  F,  M,,  +  F,  M,i.  -  Fi.  M,,]/D  (D.6b) 


S3=   [F1   M,3  +  F2M1^F3M11   +F4M12]/D  (D'6c) 

S4  =    [-   F,    H]k  +  F2  M,3    -   F3   M,2  +  ^  M, ,  ]/D.  (D.6d) 


121 


Since  there  are  only  four  unique  constants,  all  the  minors  found  in  (D.2) 
can  be  expressed  as  one  of  four  possible  minors.   The  solutions  to  the 
four  minors  present  in  (D.6)  are 

Mn  =  C,  (C,2  +  C22  +  C32  -  Ck2)    +  2  C2  C3  C^  (D.7a) 

M]2  =  -  C2  (C,2  +  C22  -  C32  +  C^2)  -  2  C,  C3  Ck  (0.7b) 

M)3  =  C3  (C,2  -  C22  +  C32  +  Ck2)    +  2  C,  C2  C^  (D.7c) 

M,i,  =  "  Ck   ("  Cl2  +  C22  +  C32  +  C^2)  "  2  Cl  C2  C3'      ^D-7d) 

The  second  important  result  of  this  work  is  the  fact  that  the  general 
solutions,  (D.6)  and  (0.7),  also  hold  for  both  the  second-order  and  third- 
order  systems  of  equations.   Naturally,  the  values  of  C.  -  C.  and  F.  -  F. 
are  different  for  each  order  of  equations. 

For  the  second-order  system  of  equations  the  solutions  S,,  S2,  S^ , 
and  Sk  are  identified  with  the  second-order  constants  I,  J,  L,  and  M 
respectively.   The  other  terms  are 

C,  =  sin(BQ2,  Z  )  sin(2  BE,  Z  )  (D.8a) 

1         b  a         D  a 

C2  =  cos(B0_2b  Za)  sin(2  BEb  Za)  (D.8b) 

C3  =  sin(BQ2b  Za)  cos(2  BEb  Za)  (D.8c) 

C.  =  cos(B0_2L  Z  )  cos(2  BE.Z  ) ,  (D.8d) 

H  b  a  b  a 
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while   the   forcing    functions   are  given   by 


Fl    =    £Yla    (Za)    cos(2  xa   Za>    "  Y2a    (Za>    sin(2  xa  2a)]/(2   BEb) 


"  Y,b   (Za}    C°S(2   BEb   Za}    "  Y2b   (Za}    Sin(*   BEb   Za>  (D'9a) 


F2  "    [Y3b    <Za>    sin<2   BEb  Za>    +  \b    (Za>    COS^2    BEb  Za>  ]/B(>2b 
+   (a/k)2    (Va/Vb)     [Y3a    (Za)    sin(2   xa   1Q)    +  Y4a    (Za) 

•    cos(2  xa   Za)]/BQ2b  (D.9b) 


F3  -    [Y|a    (Za)    sin(2   xa  Zg)    *  Y^    (Za)    cos(2  xg   Za)]/(2   BEb) 

"  V,u   (ZJ    sin(2   BE.    Zj    +  Y-h   (Zj    cos(2   BE.    Zj  (D.9c) 


lb   ^a> 


[Y3b   (Za)    cos(2   BEb  Za)    -  Y^    (Za)    sin(2   BEb  Z^/BQ^ 


(«A)Z    (Va/Vb)     [Y3a    (Za)    cos(2   xa   Za)    -   Y4a    (Za) 


•    sin(2  xa  Z   )]/BQ2b. 


(D.9d) 


For  the  third-order  solutions  the  constants  Cj  -  C^  are  identical 
with  the  first-order  values  given  by  (D.3).   The  solutions  S],  S2,  S^ , 
and  S^  are  identified  with  the  third-order  constants  P  ,  Q  ,  R.,  and 
S.  respectively.   The  third-order  forcing  functions  are 
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Fl  ■  -  Y5b  (Za)  C°s(BEb  Za)  "  Y6b  (Za>  sin(BEb  Za)  +  [Y5a  (Za> 

•  cos(*a  Za)  +  Y6a  (Za)  sFn<xa  Za)]/BEb  (D-,0a> 


F2  =  (Va/Vb)  [Y7a  (Za)  sin(xa  Za)  +  Yga  (Za)  cos(xa  Za)]/BQ1b 

"  [Y7b  (Za}  Sin(BEb  Za}  +  Y8b  (Za}  C°s(BEb  Za)]/BQ,b   <D-10b) 


F3  =  [Y5a  (Za>  S'n(xa  Za>  "  Y6a  (Za}  C°s(Xa  Za} ]/BEb 

"  Y5b  <Za>  sin<BEb  Za>  +  Y6b  <Za>  COS<BEb  Za>  (D-,0c> 


F4  -  -  (Va/Vb)  [Y?a  (Za)  cos(xa  Za)  -  Yga  (Za)  sin(xa  Za)]/BQlb 

+  [Y?b  (Za)  cos(BEb  Za)  -  Ygb  (Za)  sin(BEb  Z^J/BQl,,.   (D.lOd) 
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